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The problem of kinematic effects in gluon and color dipole cascades is addressed in the large N c 
I ^ ■ limit of SU(N C ) Yang-Mills theory. We investigate the tree level multi-gluon components of the 

gluon light cone wave functions in the light cone gauge keeping the exact kinematics of the gluon 
emissions. We focus on the components with all helicities identical to the helicity of the incoming 
gluon. The recurrence relations for the gluon wave functions are derived. In the case when the 
virtuality of the incoming gluon is neglected the exact form of the multi-gluon wave function is 
obtained. Furthermore, we propose an approximate scheme to treat the kinematic effects in the color 
dipole evolution kernel. The new kernel entangles longitudinal and transverse degrees of freedom 
I 1 and leads to a reduced diffusion in the impact parameter. When evaluated in the next-to-leading 

Oh' logarithmic (NLL) accuracy, the kernel reproduces the correct form of the double logarithmic terms 

of the dipole size ratios present in the exact NLL dipole kernel. Finally, we analyze the scattering of 
the incoming gluon light cone components off a gluon target and the fragmentation of the scattered 
state into the final state. The equivalence of the resulting amplitudes and the maximally-helicity- 
I , violating amplitudes is demonstrated in the special case when the target gluon is far in rapidity 

■ from the evolved gluon wave function. 

on : 

ON ■ I. INTRODUCTION 

^ ; 

In this paper wc investigate the possibility of using information available in the exact tree level QCD amplitudes in 
, the large N c limit to improve the evolution equations of QCD. Following the Mueller's color dipole model approach 
P, 0] we choose the framework of the light cone perturbation theory (LCPT) d, d, [j| Q and adopt the light cone 
r". \ gauge. The attractive feature of the color dipole model is its probabilistic interpretation in terms of the multi-dipolc 
densities. This allows to extract the total cross-section for a scattering of two small sized onia (quark-antiquark 
systems) at a large energy, -y/s, from the tree level amplitudes in the soft-gluon limit. The color dipole model provides 
a resummation scheme of the leading logarithmic (LL) contributions in powers of a s logs. It is remarkable, that the 
analysis of the tree level amplitudes in the color dipole framework reproduces the results of the leading logarithmic 
Balitsky-Fadin-Kuraev-Lipatov (BFKL) evolution equation Q, that incorporates quantum loops. This suggests that 
a more accurate treatment of the color dipole (or gluon) cascade at the tree level but beyond the soft gluon limit, 
may provide a lot of insight into higher order corrections to the BFKL equation. 

The basic object in our analysis is the light cone wave function of a gluon, that may be decomposed into multi-gluon 
components. The whole analysis of this paper is carried out in the large N c limit. As the simplest starting point 
we choose the configuration in which both the initial state gluon and all the gluons in the wave function carry the 
positive helicity We aim at obtaining an exact description of the gluon cascade, so we keep the complete information 
about the gluon kinematic variables along the cascade. This implies that both the form of gluon splitting vertex and 
the energy denominators of the light cone perturbation theory are treated exactly. This is different from the soft limit 
where the approximate forms were used. 

As the first result, we find an exact form of a tree level multi-gluon components of the real gluon wave function in 
the chosen helicity sector. This is achieved by solving recurrence relations for the multi-gluon components of the gluon 
light cone wave function in the momentum space. The natural variables in this formulation are closely related to the 
spinor products which are used in the construction of the maximally-hclicity- violating amplitudes (MHV), see e.g. 
[H, Q . This wave function can be also analyzed in the transverse coordinate space, after a suitable Fourier transform 
is performed. Curiously enough, the obtained wave function in the coordinate space coincides with the one obtained 
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in the soft gluon limit. This coincidence, however, is expected to hold only for the particular choice of hclicitics, that 
has been assumed. 

In general, the description of the gluon cascade within LCPT is non-local, the energy denominators involve momenta 
of all gluons, and that makes the resummation difficult in a general case. In order to obtain a more practical 
improvement of the dipole model kernel, we propose a local approximation to the 1 — *■ 2 gluon splitting kernel in 
which the value of energy denominator is obtained from the invariant mass of the produced gluon pair. This leads to 

a new color dipole kernel. The form of the kernel in the soft gluon limit ^f- d 2 x 2 Ji a ^i gets modified to 
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where x 01 is the parent dipole vector, x 02 , x 21 , are the daughter dipole vectors, z is the fraction of the longitudinal 
momentum of the softer daughter gluon, and K\ is the modified Bcsscl function. The modified kernel becomes 
equivalent to the leading logarithmic kernel for zx\ 2 "C ^oi an( ^ zx \\ ^ x \\ but differs significantly otherwise. In 
more detail, the production of larger dipole sizes is exponentially suppressed above the cutoff size which depends on 
the longitudinal momentum fraction of the soft gluon. The modified kernel introduces corrections at all orders beyond 
the leading logarithmic approximation. In particular, we show that the modified kernel reproduces an important part 
of the NLL kernel, namely the double logarithmic terms. Interestingly enough, these are the dominant terms that 
violate the two dimensional conformal symmetry of the BFKL-dipole kernel both in N = 4 Supcrsymmetric Yang- 
Mills (SYM) theory and in QCD at NLL approximation [IS El El EE G3 ■ The suppression of the large dipole 
emission implies also a suppression of the diffusion in the impact parameter. 

Finally, we consider the production amplitude of n gluons in a scattering process of two gluons. We impose additional 
conditions that the produced state is arbitrarily far away in rapidity from the target gluon, and that both the incoming 
and all the outgoing gluons carry positive helicities. Using the LCPT we reproduce formally the Parke- Taylor form 
of the MHV amplitudes Q in this limit. In the intermediate steps of the proof we analyze in detail the fragmentation 
amplitude of m virtual gluons into n > m real gluons. We find that an apparent entanglement of the fragmentation 
process of different gluons in LCPT can be unwound at the tree level, and in fact, fragmentation amplitude of each 
of the gluons is independent. We find an exact form of the tree level amplitude for a gluon fragmentation, again 
assuming that the hclicity of the incoming and outgoing gluons arc positive. We observe, that the n-gluon component 
of the incoming real gluon wave function is dual to the fragmentation amplitude of a gluon into n real gluons i.e. the 
latter has the same form as the former under replacement of the gluon transverse positions by the gluon transverse 
(light cone) velocities. 

The structure of this paper is as follows: in the next section we present the overview of the color dipole model and 
the next-to-leading order corrections to BFKL. In Scc. lIIII wc set up a general formalism for computing the multi-gluon 
tree-level wave function with exact kinematics in the light cone perturbation theory and in the light cone gauge. We 
derive the recurrence relations for the wave functions with different number of external gluons in momentum space 
and in transverse coordinate space. In Sec. IIV1 we derive the modified kernel for the dipole evolution in transverse 
coordinate space. We show that when expanded to the next-to-leading order in the strong coupling constant the 
double logarithmic terms in the next-to-leading calculation are recovered. We discuss the implications for the impact 
parameter dependence of the cross section. In Sec. [V] we construct the scattering amplitudes for external on-shell 
partons and demonstrate the general factorization property of the fragmentation in the light cone perturbation theory. 
The equivalence with the MHV amplitudes is also demonstrated in the high energy limit. In the last section we briefly 
summarize our results. The Appendix contains the most technical part of the calculations. 



II. OVERVIEW OF COLOR DIPOLE EVOLUTION IN QCD 



In a conventional approach to the high energy limit one computes the amplitude for scattering of two highly energetic 
particles which interact via exchange of gluons in the ^-channel. The emitted gluons in this scattering process are then 
resummed in the multi-Rcggc kinematics. The amplitude turns out to be governed by the exchange of a Reggeon with 
the quantum numbers equivalent to that of the vacuum: the Pomeron. Its intercept is evaluated from the solution 
to the Balitsky-Fadin-Kuraev-Lipatov (BFKL) Q evolution equation at small x. This equation describes the change 
of the amplitude with increasing rapidity between the scattering initial particles. In the groundbreaking paper [l| 
Mueller demonstrated that the small- x cascade of gluons can be alternatively described as the evolution of the hadron 
wave function. The subsequent emissions of the soft gluons in the hadron wave function can be computed in the light 
cone perturbation theory. By soft gluons we mean here gluons whose longitudinal momentum fractions are much 
smaller than that of the incoming particles, the transverse momenta of the gluons are not restricted. In the large N c 
limit it was demonstrated that the correct degrees of freedom at high energy are qq dipoles. Furthermore, by taking 
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the Fourier transform with respect to the transverse momenta it was shown that the soft gluon emissions factorize in 
the transverse coordinate space. The resulting evolution equation for the scattering amplitude was then formulated 
by considering the interaction of the hadron wave function with the target. It can be written in the following form 

dNy(x 01 ) _ ^ fd\ 02 x} n [Ny{M ) + Nvix^) - N Y { M )] . (2) 



dY J 2n xf 2 ^ 



-U2 



Here, Ny is the scattering amplitude for a dipolc on the target, Y is the rapidity for the process and x. 0: x 1: x 2 are 
the two-dimensional vectors which specify the position of the color dipoles. To obtain the solution to this equation 
one needs to specify an amplitude at the initial rapidity Ny . The solution to this equation for the dipole scattering 
amplitude turned out to be equivalent to the one from the original BFKL equation, at least in the case of the inclusive 
quantities [H, like the total cross section. 

The dipole evolution proved to be a powerful framework for investigating the unitarity corrections [lH fioj ] in the 
scattering processes. By taking into account multiple rescattcrings of the gluon components of the evolved hadron 
wave function onto the nuclei the non-linear evolution equation was derived [l7l . [l8l ] . It differs from ([2]) by the presence 
of an additional nonlinear term 

-Ny^x^Ny^x^) . 

which tames the growth of the amplitude with the energy and leads to the gluon saturation. The same equation was 
derived in the context of the high energy operator expansion with Wilson lines [l9|, [2(| HH, H3, HI, H3, Hll and in the 
Color Glass Condensate framework [Iff, [27], HH, see also [2j| for a nice review. It is called the Balitsky-Kovchegov 
equation. It is important to note that the kernels both in the linear and non-linear equations are exactly the same 
(to the NLL level in the dipole approach). 



A. Next-to-leading corrections and the kinematical effects 

Similarly to the t-channel approach the dipole evolution kernel receives large corrections at next-to-leading order. 
The next-to-leading order kernel in the momentum space was first evaluated in [3(il l3ll . l32l . [H[ in the forward case, 
and later in [34], l3a ] in the non-forward case. The dipole kernel at next-to-leading accuracy was then computed in 
[ToL l3l| via direct calculation and also in [ll|, [T^, 0, 13 by performing the Fourier transform of the existing result for 
the non-forward BFKL kernel in the momentum space. The running coupling corrections to the dipole kernel were 
also evaluated in [37|]. There are different sources of the next-to- leading corrections and it is quite well known that 
the major part of them which is common to QCD and INT = 4 Super Yang Mills theory comes from the kinematics, 
see for example [H, [H, E3- The slow convergence of the small x series can be argued intuitively as follows. In the 
standard perturbative collinear approach one identifies the large scale /i 2 which characterizes the short distance part of 
the process with strong interactions. The limit of large /i 2 implies automatically the smallncss of the strong coupling 
constant, which allows to seek the solution to the hard scattering cross section and to the splitting function in the form 
of the perturbative series. This series should be convergent at least in the asymptotic sense. In the high energy limit 
the large parameter is the total available energy y^s, the strong coupling however is not naturally small in this limit, 
and can take large values too. Therefore there is no guarantee that the expansion in powers of the strong coupling 
in the high energy limit will be st rong ly convergent. Indeed, it is observed that the small x expansion is poorly 
convergent, see for example [40L l4ll. [42| . By assuming the multi-Regge kinematics one makes strong approximations 
onto the phase space of the produced particles. As a result the perturbative methods which rely on the expansion in 
terms of the coupling are simply inefficacious in correcting these kinematic approximations at higher orders. 

Even before the next-to-leading correction to the BFKL kernel was computed, it was shown that the constraints 
from the more careful treatment of the kinematics (while being formally of the higher order) are numerically very 

Taking this effect into account leads 



important and significantly reduce the growth of the gluon density [38. _ 

to the predictions which arc in a good agreement with phenomenology [43l. 44|, [4f| . It was first demonstrated by 
Salam [4ll . l4r| that the kinematic effects manifest themselves at next-to-leading order in the form of the double 
logarithms of ratios of transverse momenta, and correspond to the higher order poles in the Mellin space co njug ated 
to the transverse momenta. The formalism of the collinear resummation developed later [4?], [H, [49|, [H(J HH, l52| was 
shown to stabilize the small x series by imposing onto the kernel the kinematic constraints and the momentum sum 
rule. The small x expansion is also stabilized by the appropriate implementation of the running coupling effects, in 
the case of QCD. Another approach to match the DGLAP and BFKL resummations was constructed in parallel in 
Refs. Pl^,|5lj5fljp|58t and similar results were found. Finally, the same problem was addressed in yet another 



scheme in Ref. 59^ 6(f 



It is therefore quite urgent to address the problem of the better treatment of the kinematics in the context of 
the dipole evolution. As the next-to-lcading corrections take into account only part of the kinematic corrections it 
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is necessary to investigate the origin and effect of the kinematical constraints on the dipole evolution. Such effects 
necessarily go beyond the fixed order calculation and are part of the all-order (in powers of a" In'™ s) resummation. 

III. EXACT KINEMATICS IN THE MULTI-GLUON WAVE FUNCTION 

A. Setup of the problem 

In this section we shall derive the exact form of the n-gluon components of the real gluon wave function, n = 
1,2,..., using the LCPT for a special choice of helicitics. Thus, we shall assume that the incoming gluon has the 
positive helicity and analyze the case when all the gluons in the wave function have all positive helicities as well. Such 
choice of the helicity state is dictated by an expected simplicity of the result. This choice also prepares the ground 
for the comparison with the Parke- Taylor amplitudes that we shall make in Scc.fV] Wc understand the computation 
presented in this section as a demonstration of a potentially useful calculational technique and the first step towards 
the exact determination of the gluon wave function with arbitrary helicitics at the tree level. 




FIG. 1: The multi-gluon wave function. The vertical dashed lines symbolize different intermediate states where we need to 
evaluate the energy denominators. It is understood that the wave function scatters finally on some target. 

To be specific, we start with a gluon with four momentum P, and the color index a, that develops a virtual fluctuation 
into states containing n gluons with momenta (fci, . . . k n ) and color indices (oi, 02, ... , a„), correspondingly see Fig. [T] 
For the initial virtual gluon with virtuality — Q 2 and a vanishing transverse momentum one has in the light cone 
variables. 



2P+ 



,Q) , 



with P ± = -^(P ± P 3 ). After n - 1 splittings the wave function is shown in Fig. [T] The last n gluons have the 

corresponding momenta labeled ki, . . . , k n , as in Fig [I] Each of this momenta has components fcj 1 = (ziP + , k~ , fcj, 
with Zi being the fraction of the initial P + momentum which is carried by the gluon labeled by i and k i being the 
transverse component of the gluon momentum. The rules of the light cone perturbation theory [H, 0, 0, require 
to evaluate the energy denominators for each of the intermediate states in Fig. [1] The energy denominator for the 
n gluons in the intermediate state is defined to be a difference between the light cone energies of the initial and the 
intermediate state in question. For the wave function shown in Fig. [T] we assume the last intermediate state is with 
?i gluons as depicted. The corresponding energy denominator for this state reads 
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where we have used 

1 2 Zl P+ ' 

and we introduced the auxiliary notation for the (rescaled) denominator D n . 

In what follows, we shall focus on the color ordered multi-gluon amplitudes, that is the amplitudes decomposed in 
the basis of color tensors, T" 1 " 2 -"" = tr (t ai t a2 . . . t a ") , where 01, a 2 , ■ ■ ■ , a n are the color indices of the gluons. Up 
to the leading terms in the large N c approximation, these tensors are orthogonal, see e.g. Q for a detailed discussion. 
This decomposition has a transparent graphical representation in 't Hooft's double line notation (see e.g. [u"l| for a 
review). The n-gluon amplitudes in the basis T aia2 --- a ™ are then represented by planar diagrams built of the effective 
colored fermion lines. This, in turn, leads to the color dipole representation proposed by Mueller [l[, in which the color 
coefficient in the triple gluon vertex is equal to unity. In what follows, we shall use the n-gluon components of the 
incoming particle light cone wave function ^ n (ki, k2, ■ ■ • , k n ), assuming implicitly, that the corresponding color tensor 
is given by the planar diagram in the double line notation, that corresponds to the definite ordering (ai, a 2 , ■ ■ ■ , o, n ) 
of the gluon color indices. 



B. Recursion relations 



We focus on components of the wave function in which all the gluons have positive hclicitics. We work in the 
light-cone gauge, 77 • A = 0, with vector n = (0, 1,0) in the light-cone coordinates. This choice of gauge defines the 
polarization four-vectors of the gluon with four-momentum k 

e(±) = e < ±) + s ^_i, i (4) 

where = (0, 0, eW), and the transverse vector is defined by = =F^(l,±i). The explicit projections of the 
triple and quartic gluon vertices show on this helicity states show that @: 

1. The helicity transitions ( ► H — h) and (H ► ) arc forbidden for the triple gluon vertex. 

2. The 1 — * 3 gluon transition given by the four-gluon vertex leads to mixed helicity composition of the three 
outgoing gluons. 

These properties indicate, that the state after splitting always contains at least one gluon with the helicity of the 
parent. Thus, if we require the final state to be composed from gluons with positive hclicitics, there can be no gluons 
with negative helicities at any intermediate step. So, the four gluon- vertex does not contribute to the splitting process 

that we consider. In fact, the whole branching process is driven by a (H ► H — h) helicity projection of the triple gluon 

vertex. 

We define the projection of the triple gluon vertex on the (±) polarization states, assuming that all the gluons are 
outgoing, 

k 2 , k 3 ) = e< Al W^W^V£££(*i> **. fe) > (5) 

where a±, 0,2, 03 are of course color indices. As we already stated, vertex V^ 1 "^ (k\, k 2 , k 3 ) vanishes if all the outgoing 
gluon helicities are the same, Ai = A2 = A3. When one helicity, say Ai = —1, is different than the others A2 = A3 = +1, 
the vertex takes the following form in the light-cone variables, 

V^ aa {ki,k 2 ,k 3 )=9f aia " 13 ■ - ^) , (6) 

\Z2 Z 3 J 

where the S- functions related to the conservation of the transverse and longitudinal '+' components of the momentum 
are implicit. Here g is the strong coupling constant and f a ^ a ^ a a ; s structure constant for the SU(N C ) color group. 

For the case of interest, of (H ► ++) transition, the amplitude is described by ([6]), with z\ being the fraction of 

the + component of the momentum of the incoming gluon. Note that all dependence of the vertex on the transverse 
momenta of the daughter gluons is completely absorbed into a variable 
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This variable is a relative transverse light cone velocity of the two gluons. Interestingly enough, the same variable is 
present when we consider the change of the energy denominator due to the splitting. In a general situation, when the 
gluon with momentum k\ belongs to a virtual gluon cascade, the energy denominator before the splitting of gluon 1 
can be written as 

k 2 

D n = D n/1 + =i , (8) 
z\ 

k 2 

where D n n = Q 2 + J2i>i T~ does not contain the energy of gluon 1. Note that, we are using here the definition Q 
for the energy denominator, which is different by the sign and with the P + dependence factored out. After the gluon 
splits into two gluons with momenta ki and k 3 we have 

7 2 ,2 

D n+1 = D n/1 + =2 + =3 . (9) 

Z2 Z 3 

In the light cone perturbation theory the transverse and the + components of the longitudinal momenta are conserved 
in the vertices therefore we have that that k l = k 2 + k 3 and z\ — z 2 + z 3 . Using this fact one can express the change 
of the energy denominator as, 

Z2Z3 ( k 2 k 3 



D n+1 - D n = — - — ■ (10) 

Z2 + Z 3 \ Z 2 Z 3 J 

It is convenient to introduce a variable that depends on the longitudinal degrees of freedom 

6s - • (ID 

Z2 + Z 3 

We therefore see that the change in the denominator ([TOjl due to the splitting is expressed through the variable i>23 
as well 

D n+1 -D n =6 3 w| 3 . (12) 

In the light-cone formulation of the QCD, the intermediate line that carries longitudinal momentum fraction Zj is 
multiplied by \j ^fz[. It is therefore convenient to follow the convention by @ and absorb such factors into all the 
gluon lines incoming and outgoing from the triple gluon vertex. Thus, we shall use 

KTx!(kuk 2 ,k 3 ) = -^V^( ku k 2 ,k 3 )=gr^ . (13) 

V ZIZ2Z3 VS23 

Following the discussion of the color structure of the amplitudes in Sec. IIII Al in what follows we drop the structure 
constants f^ 11 ^. 

Now, we can collect the vertex and the energy denominator together to get the effect of the gluon splitting on the 
virtual multi-gluon state wave function. We shall introduce the following notation. Let \l/ n _i(fci, k 2 , ■ ■ ■ , fcjj+ii ■ • ■ , k n ) 
(with fcji+i = ki + fcj+i) be the (n — l)-gluon wave function in momentum space before the splitting of gluon with 
momentum feft+i, and ^ n (ki, k 2 , ■ • • ,ki, fci+i, . . . , k n ) the wave function after splitting of this gluon. Then the splitting 
of the gluon with momentum k ii+ i acts as follows, 

ty n (k\, . . . , fcj, fcj+i, . • . , k n ) = — j= — - " +1 2 if? n —\(ki, . . . , fejj_)_i, . . . , k n ) , (14) 

where the color degrees of freedom are treated in the way described in Sec. IIII Al 

The formula Q14p is the recurrence prescription for obtaining the wave function with n virtual gluons from the 
wave function with n — 1 gluons. Of course to obtain the full recurrence formula one needs to sum over the different 
possibilities of the splittings which will give us the following result 

n (_) 

*n(/fl,/%2) •••,&*») = "FT ^ ~ n^-= ~ ^n~l( fc l: ■ ■ • M-lf, ■ ■ ■ ,k n ) , (15) 
Un i=2 V«"li 



where D n is the denominator for the last intermediate state with n gluons. The form of the splitting given by |[T7 
and P3)l looks highly symmetric, and all dependence on momenta of daughter gluons i and z — 1 is embedded into two 
variables: j and v_i_i j. In this calculation we have kept the exact kinematics and therefore the energy denominator 
includes the full dependence on the momenta in the whole cascade. This is understandable, as in the case of the exact 
kinematics the new wave function has to carry the full information about the old wave function. We will come back 
to this point later in Sec. Ill when we discuss the high energy limit. 
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C. Relation to helicity amplitudes and the collinear limit in the on-shell case 

It turns out that the variables Vj k that we used to construct the wave functions in the previous subsection are 
related to the variables used in the framework of helicity amplitudes, see @ for a nice review. Namely, for given pair 
of on-shell momenta ki and kj we have that 



where the symbols [ij], (ij) are the spinor products defined by 

(i\j) = (i-\j+), [ij] = (i+\j-) ■ (17) 
The chiral projections of the spinors for masslcss particles arc defined as 

\i±) = V±(fci) = \{l±ls)i>{h) , (±i\ = i^(k~j, (18) 

for a given momentum ki . The spinor products are complex square roots of the total energy mass squared for the 
pair of gluons 

(ij)[ij) = (k i + k j ) 2 , (19) 
and they also satisfy (ij) = [ij]*- Using the above definitions ()16|19|) we have that 



ij)[ij] = 



Zi Zt 



which is real and positive for the on-shcll gluon momenta. Finally, combining ([7]) and (|16[) we obtain 

( ij > = ^/^] e (+) • Uij , [ij] = V^j ' ^ij > ( 20 ) 

and the dependence on the transverse momenta in the light cone wave function can be expressed by ( ij ) and [ij]. 

Using these expressions we can check the collinear limit for the on-shell case. The Eq. (|T4l) is part of the recursion 
relation (TTSJ) for the off-shell multi-gluon wave function. It actually describes the situation in which the gluon with 
momentum ki 2 splits into two daughter gluons, with momenta k± and ki respectively It is interesting to investigate 
the collinear limit of the on-shell amplitude, which should get factorized. To get the on-shell amplitude one needs to 
drop the non-local denominator -D n _i in (fT4]) . The factorizable limit for gluons 1 and 2 is then 



*£ ll2) (fcl, **,...,*„) = ^ ^ n -l(fc 1 + fc 2i fc 3 ,...,fe„) = 

1 Q \12] 



/Cl2 sjz(\ - z) S12 



y n - 1 (ki + k 2 ,k 3 ,...,k n ) (21) 



where z = z\/{z\ + z 2 ) and s i2 = (ki + k 2 ) 2 ■ Relation (|2"Tj) is, modulo 1/V?i2 coefficient and the sign convention, 
exactly the factorization relation on the collinear poles for the kinematical parts of the dual amplitudes as shown in 
. It is interesting to note that the only thing that we have done here is to identify the gluons with momenta k\ , k 2 
as originating from the splitting of the gluon with momentum k\ 2 and therefore we have selected only one splitting 
out of n — 2 possible combinations. 



D. Multi-gluon wave function in the coordinate representation 

It is interesting to investigate the form of the recursion relation (fTi)) in the transverse coordinate representation. 
By this we mean performing the Fourier transform with respect to the transverse components of the momenta, just 
like in the original dipole approach [l[ . In this section we will consider a special case of the wave function where in 
the initial state we have only one gluon. In the original approach the initial state was the quark-antiquark pair. We 
will discuss the latter case in the next section where we will derive the modified dipole kernel. In this section we want 
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to demonstrate calculational techniques which allow us to resum the consecutive splittings in the wave function. We 
therefore assume that initially we have the just one gluon, which has momentum k\. The initial wave function in the 
transverse space is then written as follows 



$!(!) = $ 1 (z 1) r 1 ) = J ^e^M,^,^) 



(22) 



where is the coordinate in the transverse space. In the rest of the paper the symbol $ will always denote the wave 
function in the momentum space, whereas $ will denote the wave function in the coordinate space. For the wave 
function with n gluons we define 



^ n (zi,r 1 ;z 2 ,r 2 ; . . -,z n ,r n ) 

d 2 k x d 2 k 2 d 2 k n 
(2tt) 2 (2tt) 2 ' ' ' (2^)2 



exp(ifc 1 ■r 1 +ik 2 -r 2 + ... + ik n ■ r n ) * n (/ci, . . . , k n ), (23) 



and similarly for and \& n _i. For the purpose of the subsequent calculation it will be convenient to change 

the variables for the transverse momenta from (ki,k 2 ) to K_ 12 = k_ x + k 2 and k 12 = £12^121 and the longitudinal 
variables from Zi,z 2 to z\ 2 = z\ + z 2 , and £12. These are actually the c.m.s. variables for gluons 1 and 2. In 
addition one has: d 2 k 1 d 2 k 2 = d 2 I£ 12 d 2 n 12 , and z\z 2 = £% 2 zi 2 . The scalar products can be rewritten as follows, 
kiLi + A-2^2 = Ml.yiR.yi + £12^12 1 where R 12 = 2l ^ + ^~ 2 , and r 12 = r 2 — r x . Next, we represent the energy denominator, 



1 



112 



A' 



12 , ^3 

Z12 z 3 



as 



1 r 00 

— = / dr exp(-rD„) 

i^n Jo 



(24) 



where recall that Q 2 is the virtuality of the incoming particle. The above formula for the denominator (f2"4"| may be 
regarded as the Hamiltonian representation with the Euclidean time r. Using this representation we arrive at 



$„(l,2,...,n) 



d 2 Kin d 2 K, 



(2tt) 2 (2tt) 2 ' ' ' (2tt) 



d 2 k n ., . )u 12 

cxp(zK 12 ■ I 12 + tK 12 ■ i< 12 + ... + « 



12 



x 



/>oo 

/ dre- rDn V n -i(K 12 ,k 3 ,...,k n 
Jo 



— —ri2 



d 2 K 12 d 2 K 12 



£l2^ J (27T) 2 (2^) 2 (2tt) 



d 2 fc 

■ ■ 7^ cxp(iK 12 • r 12 + iK 12 -R 12 + ... + ik n ■ rj 



dr exp 



112 
£12 



-2 
-12 



Z12 



/ dM -^- 1 ex P H fe + fc2)-ri-^-Zl2 



ife n -d-i)^-i(l',- 1)0 , (25) 



where we have used the recursion formula (|14[) and the inverse Fourier transform for the n gluon wave function. Note, 



that the primed z values are the following: z[ = Z\ + z 2 , z' 2 = 23, 
momenta are Gaussian and can be easily performed, which results in 



z n . The integrals over all transverse 



*n(l,-- 
x exp 



,n = % 1 



£l2V?12 



62^12-23 



dr 







f— y 

V 47TT / 



.d 2 r> 



[62£ 2 2 + z 12 (r! 1 - R l2 f + z 3 (r! 2 - r 3 ) 2 + . . . + z n (r[ n _ x - r n ) 2 ] - Q 2 r 



$„_!(!', 2',..., (n- 1)') 



(26) 



The integral over the 'time' r can be performed to give the expression with the modified Bessel function Ki 



*„(!,..., n) 



f(-) ■ 



^12 



Z\Z2 . 



dV 1 ...d%_ 1 $„_i(l',2' 



(n-1)') 



(27) 
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where 

A = a 2 r 12 + z 12 (r' 1 - R 12 ) 2 + z 3 (r£ - r 3 ) 2 + . . . + ^(r^ - rj 2 . (28) 

The integral kernel (— ^^-)" -Ko f \/ Q 2 ^J can be rewritten into a more elegant form using the relations between 
the modified Bessel functions and their derivatives, see for example [62j 

[x- n K n {x)\ = (-l) m x- n ~ m K n+m (x) , 

which gives 

(-;£)" -p^(t) , *"<^ 3 )- 

Using the above relations we can recast the recurrence relation (|27|l into 

*n(l, •••,») = i ^/^ ZlZ2 --- Zn j d ~^(0^ (t) 2 *«(V^) 2 '-- ■•. («-!)')> (29) 

with A defined above (|28|) . The formula (|27|) is the prescription for the off-shell tree level wave function with exact 
kinematics. Since it depends on the coordinates of all n — 1 gluons through the variable A defined above, it is quite 
complicated. Also, it should be kept in mind that this is just formula for one particular splitting, one needs to sum over 
all possible splittings like in ([15]). We will show nevertheless that in the case where the incoming particle is on-shell 
there are significant simplifications, which allow to resum the multiple gluon splittings. The crucial difference with 
respect to the leading logarithmic approximation with the Regge kinematics is the appearance of the modified Bessel 
functions K n which contain the information about the gluon splitting. In the original approach the splitting of the 
gluon leads to the expression which is just a power in the transverse coordinates. This translates into the power-like 
behavior of the splitting kernel in the dipolc equation. Here, because the kinematics is kept exact the functional 
dependence is governed by the Bessel functions, which for large values of their arguments have exponential behavior 
asymptotically. This will result in a qualitative difference when investigating the impact parameter dependence of the 
scattering amplitude. We will come back to this problem and discuss it in more detail at the end of Sec. IIVI 




E. Resumming the multi-gluon wave function in the case of the on-shell incoming gluon 

We consider here the multi-gluon wave function that originates from subsequent splittings of an on-shell incoming 
gluon with helicity +. One can also alternatively think about it as the incoming particle with a large momentum 

P + such that P~ = —^p+ is very small, at least as compared with the particles in the wave function. This will 
result in energy denominators which do not contain the initial P~. We will assume that all gluons have + helicities 
which should be a situation in the high energy limit, where the helicity flips are suppressed. We shall use the complex 
representation of the transverse vectors: Vij = e> + > ■ v^, v*j = (f~> ■ v^, and a useful notation, 



J (lH2-. ip)(jl]i...jq) 



fcij + ki 2 + . . . + ki p kj 1 + kj 2 + . . . + kj q 

Z%\ ~t~ ^22 ~t~ ■ • ■ 4" Zi p Zj l -}- Zj 2 + Zj q 



(30) 



_ (Zn + z l2 + . . . + z ip )( Zjl +z 32 +... + z jq ) 
P Zn + z i2 + . . . + z. lp + z h + z J2 + ... + z jq 

with notation ki = e' + ' • k i . The global momentum conservation (5-functions, A'™) = 5^ ( Y^i=i —i ) ^ ( 1 ~ z i )i 

that are present in all the expressions for the wave functions will not be displayed explicitly. Thus, the incoming state 
has the wave function, 

*i(l) = l. (32) 

In the following discussion we will consider color ordering in the amplitudes, therefore we will suppress color degrees 
of freedom. In general, for a color-ordered amplitude, the gluon splitting acts on the wave function as derived in 
Eq. (fT51) (for on-shell initial state though) 

-D n+ i* n+ i(l,2,...,n + l) = 



10 



J nn+1 



= ff ^p=* n (12,3,...,n+l) + 0-^p=* n (l,23,... J n + l) + ... + g—L- 

VS12 V?23 V^nn+i 



* n (l,2,...,nn + l), (33) 



with D rl+ i = fci/^i + fc^/^s + • • • + fc^ +1 /z„+i. We have introduced the notation * n (l, . . . , i — 1 i, . . . ,n + 1) where 
i — 1 i means that it is the gluon with the momentum ki—n = + fe,-. After the first splitting one gets 



*a(1.2) = -9 



12 



1 1 



/fl2 62kl2| 2 \/£l2 C12W12 ' 

where we have taken that P~ = 0. According to (|3"3"|) . the next splitting leads from \&a(l> 2) to ^3(1, 2, 3): 



(34) 



-1)3*3(1,2,3) = -9 



:* 2 (12,3) 



"23 
^23 



* 2 (1,23) 



12 



■23 



v/£ 12 £(12)3 £(12)3 «(12)3 V^23?l(23) Cl(23) «1(23) 

This expression may be, after some simple algebra, simplified using £i2£(i2)3 = ^23^1(23) = z ^ z ^+ z 
obtains, 



* 3 (1,2,3) = g 2 



1 



1 



1 



£(12)3£l(23) w (12)3 ^1(23) 



(35) 

Z1Z2Z3. One 
(36) 



Note that, the energy denominator D 3 disappeared from the equation as it has canceled with the numerator when 
finding the common denominator for expression (|35[) . The same procedure can be iterated further. We shall give 
below the explicit form of the wave function obtained for 4 gluons and then present a generalization to an arbitrary 
n. Thus, for n = 4 we found: 



* 4 (1,2,3,4) 



1 



1 



1 



-9 



1/Z1Z2Z3Z4 £(123)4 £(12) (34) £l(234) ^(123)4 u (12)(34) ^1(234) ' 

and for a general integer n > 2 one expects, 

1 1 



(37) 



n—1 n— 1 



* n (l,2,...,n) = (-!)"" V 



V z l z 2 ■ ■ ■ Z n £(12...n-l)n £(12...ri-2)(n-l ri) ■■■ £l(2...n) 



«(12...n-l)n «(X2...n-2)(n-ln) •■• w l(2...n) 



(38) 



This formula was explicitly verified for n — 2,3,4,5,6. The proof for arbitrary n can be done by mathematical 
induction and proceeds as follows. We assume that the wave function *„ satisfies the above conjecture (|38|) . Using 
(|33| the wave function ^ n +i has then the form 



D n+1 H> n+1 (l,2,...,n+l) = gJ2 -%^=*n(l,2,...,(ii + l),... ) n+l) 

f=l V&i+l 

n » 

(_i)«-y £ _^»±i 



1 V£ii+1 V z l z 2 ■ ■ ■ (^i + Zi+l) • ■ ■ Zn (£(12...n)n+l £(12...n-l)(nn+l) ■■■ £l(2...n+l))' 

7 " v- ( 39 ) 

(V(12...n)n+1 f(12...n-l)(nn+l) ■■■ "l(2...n+l) j 

We have inserted the symbol ' to denote the fact that the indices ii + 1 have to be taken together, or in 
other words in each term of the sum in ([55)1 for a given i the denominator does not have the term of the form 
£(i...i)(i+i...»+i) u (i...i)(»+i...n+i)- The expression under the square root is 

£,ti+l ZlZ 2 . . . (Zi + Zi+l) . . . Z n = Zl * t+1 Zl z 2 . . .( Zi + Z i+ i) ...Z n = Z X Z 2 . ■ ■ Z n +1 , (40) 

Z% + Zi+i 
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because i goes from 1 to n, there arc n + 1 terms in the product (|40p . Now we write (|39[) in a form with the common 
denominator for all the terms which gives 

n \Jl (10 . i 11 r 1l»„n Z)"=l W i*+l^(l...i)(i+l...n+l) t '(l...i)(i+l...n+l) 

i'n+l ^n+H 1 ) A • • • jTl + i-) — {— J-J 5 



y/ZiZ 2 . . . Z n +i £(l2...n)n+l £(12...n-l)(« n+1) • ■ ■ Sl(2...n+1) 

• (41) 

v (12...n)n+l "(12...n-l)(nn+l) •■■ u l(2...ri+l) 

The denominator in the above equation is just what gives the wave function ^ n +i- An essential point is to show that 

n n+1 

*^2vii+l£(l...i){i+l...n+l)V(l...i){i+l...n+l) — ^ ~ = Ai+1 ' ( 42 ) 
i=l 3=1 J 

that is that the numerator on the right hand side in ([4Tj) cancels D n +i on the left hand side of this equation. We 
notice that, 



(l...i)(i+l...n)£(l...i)(i+l...n) - ( "^j— ^7 + 7Z ... + Zi )) ~ ( 43 ) 



Therefore the left hand side of (|42[) is 

n n i 71 / k k- \ 1 

y]«ii+l^(l...i)(i+l...n+l)U(l...j)(i+l...n+l) = y^ii+l]T)&j = X! T _ T^ - X] k 



1=1 3 = 1 i=l x T ' j = l 

n , i-1 



~ Zi . . .. ^2+1 ... . ~ Zi . ~ Zi . Z n +1 n Zi+i 

i=l 3 = 1 i=l 3 = 1 i=l i=2 3 = 1 2—1 3 = 1 

The second and fourth term cancel, and the first and third give 

k 2 

J2~ = D n + i- (45) 

. . Zi 

2 — 1 

which is exactly the denominator, so we have proven (|42|) . Note that, we have also used 

n 

^ ki = fcn+1 j 

2=1 

i.e. the condition which is coming from the energy-momentum conservation. This completes the proof, as ^ n +i is 
exactly of the form for n — > n + 1. 

The resummcd form (|38[) looks relatively simple. Its dependence on transverse momenta of gluons factorized into 
reciprocals of V(i2...p)(p...n) with the splitting index p takes all possible positions from p = 2 to p = n — 1. The wave 
function given by (|38[) may be also expressed in terms of the ( ij ) symbols in the following way: 

* n (l,2,...,n) = 



\J Z 1 Z 2 ■ ■ ■ Z n *\/£(12...ri-l)ri £(12...n-2)(?i-l n) ••■ £l(2...n) 



((12...n- l)n) ((12...n-2)(n-ln)) ... (1(2... n)) 



The wave function given by (|38[) can be simplified further. 
Denoting fc(i... p ) = fci + . . . + fc p we may write: 



(46) 



2, . . . ,») = (-I)"" V" 1 - — . (47) 

^/ZiZ 2 . . . Z n K(1)K(12) ■ ■ ■ K(i2...n-1) 
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One may go back to the real representation of the transverse vectors and perform the Fourier transform of (|47[) , 



d 2 k d 2 k 

$„(zi,r i; ...;z„,r„) = / ■ ■ • t^t| exp^ ■ r 1 + . . . + ik n ■ r n ) * n (zi, k 1: . . . , z n , k n ). (48) 



It is convenient to change the momentum variables: 

{^1,^2; ■ ■ -kn} — * {k(l),k(12)i ■ ■ ■ i fc(12...n)}- (49) 

The Jacobian of the transformation is trivial, and one expresses the gluon momenta in the following way: 

hl=h.(l)i k.2 = k.(12) ~ ^(1)> kn — k(12...n) ~ k(12...n-l)- (50) 

The Fourier exponent takes the form, 

exp(% ■r 1 + ... + ik n -r n ) = exp(ifc (1) • (r 1 ~r 2 ) + ik {12) ■ (r 2 - r 3 ) + . . . ik {12 „_ 1} • (r„_ x -rj +ik(i%... n ) -L n ), (51) 
and the final answer for the wave function in position space reads 

rh ( \ ( i\n-l n-lr | i \ 1 t ( ~ ) Li2t [ ~ ) L2Z ■■■^ ) L n -ln 

$„ zi,r 1 ;...;z n ,r„ = -1) g 6 1 - > = g— 5 5 • ( 52 ) 

' ^JZ\Z2 ■ ■ ■ Z n 



i=l 



Surprisingly, this is precisely the structure of the gluon wave function in the LL dipole model. However, it is crucial for 
the derivation that the incoming gluon is on-shcll, and that we chose a very special component of the wave function, 
with all positive gluon hclicities. Note that the on-shellncss (or quasi-onshcllness) of the incoming particle means that 
in the light cone perturbation theory we neglect in the energy denominators the energy of this particular particle and 
keep only the energies of the subsequent emitted gluons. 



IV. DIPOLE KERNEL WITH APPROXIMATE KINEMATICS 
A. Dipole wave-function in the leading logarithmic approximation. 

In the previous section we have presented a general approach for the treatment of the light cone wave function while 
keeping the kinematics exact. We were able to construct the recurrence relation for the off-shell case and resum the 
wave function for the special case of the on-shell incoming particle and with all the gluons having the same helicities. 
However, for the practical purposes, it would be desirable to construct the evolution equation in rapidity that would 
resum the gluon emissions just like in the original dipole approach developed by Mueller [l|. We will show in this 
section that it is possible, though not all the exact kinematical effects are kept here. The modified equation should 
include at least part of the corrections due to the exact kinematics. These corrections will prove to be sufficient to 
recover the double logarithmic corrections found in the next-todeading calculation. 

Let us therefore recall the original construction of the dipole wave function at small x with soft gluons which was 
first developed by Mueller [l[. Unlike the previous section, where we took the initial particle to be the gluon we 
will consider here the qq pair in the color singlet state, which is the onium. The heavy onium wave function, which 
consists of the quark- antiquark pair without any additional gluons, is defined as 

y {0) (zi,ki), where Zl = kf/P+. (53) 

The initial momentum is P, the quark has four-momentum k\ = (k^~ , k^ , fc x ) and the antiquark has momentum P—k\. 
In the transverse space the wave-function reads 

$ (0) (^oi) = / H^- 1 ^ 01 * (0) (-i,fci) , (54) 

where x 01 is the size of the dipole 01 in the transverse space. Strictly speaking, this representation is a mixed one: 
coordinate space in the transverse components and momentum space in the longitudinal degrees of freedom. 

Such representation was shown [l[ to be very convenient for the purpose of studying the high energy limit, as only in 
this limit there is a decoupling of the transverse and longitudinal degrees of freedom. Let us investigate the emission 
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of one gluon with momentum see Fig. [2] If we were to keep the kinematics exact we would arrive at the analog of 
the formula (TH| for the case of one splitting, where say the gluon is emitted from the upper anti-quark line (Fig. [5]) 



(i) 



gt a u{ki)^ ■ e 2 u(ki + k 2 ) 
D~i 



(55) 



Here the exact energy denominator has the form 



D 1 = P- - [(P - h 



(56) 



and similarly for the second emission from the other quark line. Following the procedure presented in the previous 
section one would obtain the recursion formulae which depend on the kinematics of the whole cascade. The simpli- 
fication that arises in the high energy limit is due to the fact that the emission of the very soft daughter gluon gets 
factorized from the rest of the wave function. In other words there is no recoil due to the emission of the gluon and 
the rest of the whole cascade is frozen and is not affected by the splitting. This allows to factorize the emissions of 
the gluons and resum their emissions in the form of the evolution equation which is differential in rapidity. We shall 
show that the leading logarithmic approach will be changed when the more exact kinematics is taken into account, 
but the resummation of the emission can be still recast in the form of the evolution equation. To this aim let us 
recall the original assumptions made in [l[ that lead to the derivation of the dipole evolution equation in the leading 
logarithmic approximation. 

The following assumptions are done in order to reproduce the leading logarithmic approximation, or the limit of 
the high energy 

• The emitted gluon is longitudinally soft: k 2 <C fcf, and one defines z\ = k^/P + ,z 2 = k 2 / P + which are the 
fractions of the longitudinal momenta. Obviously z 2 <C z\. 

• The coupling of the gluon to the quark (antiquark) is eikonal 

.gi a u(fci)7 ■ e 2 w(fci + k 2 ) = 2gt a k x ■ t 2 , 

where t a is the color matrix in the fundamental representation and e 2 is the polarization vector of the emitted 
gluon with momentum k 2 . 

• Since the gluon is longitudinally soft k 2 <C kf one can keep in graphs the leading term in the energy denominator 
(see Fig. (2) 



i 



p- 



[(P -h-k 2 



ki + k 2 ] k 2 



because k 2 is the dominant term (k 2 = ki,/2k 2 and k 2 is very small, k 2 <C k^,P + ). Here, P is the initial 
momentum, see Fig. [5J 




^1 ^2 ? — 1 — 2) 



(k 2 , k 2 



(kt,h 




(p+-k+,-k 1 ) 



(k+ -k^h-kl) 



FIG. 2: Dipole wave function with one soft gluon. Dashed line indicates the energy denominator for the intermediate state. 
If one puts together all these assumptions one arrives at the formula for the wave function with one emitted gluon 



• k 



& 1 \z 1 ,k 1 ;z 2 ,k 2 ) = 2gt a ^p-[^°\z 1 ,k 1 ) - * (0) (^i,*i + k 2 )] 9(*i - z 2 ) 



k 2 



(57) 
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At this stage there is only ordering in the longitudinal fractions z\ 3> ^2 which is extended to the weak inequality 
z i > z 2 hi the subsequent calculation. The transverse momenta are unrestricted in this approximation. One can 
perform the Fourier transform of the above expression to obtain the formula for the wave function with one soft gluon 
in the transverse space 

* (1) (*l,Soi;*,a H ) = / (2~)2 (2)2 e<fa-*»-Hfa-*» ^(^fe^fe) , (58) 
which yields the explicit result 

zVizuXo^Xn) = -^-(^■-^-)-e a ^°\z 1 ,x 01 ), (59) 

7T V^O X 21 / 

where x?j = x?j . We see the advantage of using the coordinate space representation because the soft gluons factorize 
in (|59p as demonstrated in [l[ . The modulus squared of the one gluon wave function has then explicitly the form 



<I> 



(i) 



2 (ft.jfcx) = f P- I % -A- 2 ■ (60) 



Z2 J 27T 



^02 ^12 



In this form it is particularly transparent that in the soft gluon limit the transverse and longitudinal degrees of 
freedom decouple. Therefore the transverse space coordinate representation is especially useful. The wave function 
with one soft gluon is just the wave function without any soft gluons times the branching probability. The measure 
(or branching probability) in this case reads 

; (61) 

^02 SK12 

is the dipole splitting kernel in the leading logarithmic approximation which appears in the dipole evolution equa- 
tion for the dipole scattering amplitude and it was originally derived in [l[. The same kernel appears in nonlinear 
version of the evolution equation with additional gluon rescatterings: the Balitsky-Kovchegov equation [lj], H0| • It is 
straightforward to check that this kernel is invariant with respect to the 2-dimcnsional conformal transformations. In 
deriving the equation (|57j) one uses crucial assumption about the strong ordering in the energies 

> fcf . (62) 

In the original approach [l[ one uses the assumption of the softness of the gluon with respect to the parent dipole 
and treats the transverse momenta as unordered which is equivalent to the Regge kinematics. This enables to make 
the approximations as described above and factorize the soft gluon contribution. However, for the consistency of the 

_ fc 2 _ k 2 

calculations we should keep the inequality (I62[) (where k l — -^r, k 2 = 7j£f) exact which gives 

i 2 7,2 

F>F' (63) 



with k^ > k^ ■ 

If there are more gluon emissions we will of course have 



2 

^ h+ h+ J.+ ' y > 



ki 3 kt 2 kt! 
k+ b+ b+ k 

h ii /t i 2 ft H 



and 



where the indices ii, . . . , 14 enumerate subsequent emissions along one branch of the gluon cascade. 

As a result of the ordering (f6"5| the region in the transverse momenta k^ for the gluon emission is limited 

e(^-^) + e(*i-^)e(*i-*i^r) ■ (65) 
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This means that there is a constraint which restricts the transverse momenta of the daughter gluon (labeled 2). which 
is the step function 0(k?, — Mtt)- The constraint (|63I65[) means that the momenta of the emitted gluons (those with 



k 2 ) are cutoff in the infrared. Note that, the ordering ([55]) is exactly the same as the one discussed in [63|, [6J, [65J , it 
is ordering in the fluctuation time. One has to think about the cascade as developing from the hadron side and the 
onium as a model of a hadron. 

Given the discussion above, the wave function of the onium with one soft gluon should be modified to include the 
constraint (|63p and with this modification reads 

& 1) (z 1 ,k 1 ;z 2 ,k 2 ) = 2gt a ^-f± - * (0) (zi, k x + k 2 )] Q(zi - z 2 )e(^ - ^) , (66) 

K_2 V Z 2 z l J 

where we used z\ and z 2 instead of k^ and k 2 ■ At this point it is no longer so trivial to perform the Fourier transform 
into the coordinate space since the transverse momenta are entangled now and the soft gluons do not factorize as 
before in coordinate space. In particular, there will be a shift (or distortion) of the original dipole since we encounter 
integrals of the type (if we choose to make the k x integral first in the following example ) 



J (2tt) 2 \z 2 zi 



(67) 



Therefore we see that the kinematical constraint emerges in the light-cone perturbation theory from the more 
exact treatment of the energy denominators in the graphs. In the i-channcl formulation of the BFKL Pomeron, the 
analogous consistency constraint arises, when one takes into account the fact that the virtualitics of the exchanged 
gluons are dominated by the transverse parts. This leads to the constraint on the transverse momenta of the emitted 
gluons, see for example [4(i| . 



B. Modified energy denominators in the dipole evolution 

In the previous section we have seen how the kinematical constraint emerges from the more careful treatment of 
the energy denominators. In this section we will demonstrate that one can include this effect into the dipole wave 
function. In particular we will arrive at a modified dipole evolution equation. Let us take the more exact version of 
the energy denominator which includes the energy of the parent emitter 

1 



With this modification the formula (I57p for the dipole wave function in momentum space with one gluon becomes 



fc^a, fc,) = 2jC-^[#( 2l , fc x ) - ^{z^k, + k 2 )} . (68) 

We still keep the vertex to be eikonal, the only modifications are in the energy denominator. Let us define the scale 

Q 2 = kl^=klz, z=^-, (69) 
k{ z\ 

and perform the two-dimensional Fourier transform of (|68p to the coordinate space 



Let us take the first term in the above formula 



(70) 



d 2 k x d 2 k 2 el k 1 . slol +ik 2 -x 02 ^(o) (k \ £2 ■ k 2 



£L2 

72 7 



2 9t a I ^e i ^^{k 1 ,z 1 )^-QK 1 (Qx 02 ) 



(2tt) 2 w ' l, 2^ ^ ^' x 2 

= 2gt a [ dW°)(r,z) f -^=L e *£i ' fen ~t) ^-QK X (Qx 2 ) ^2_1§02 (n) 
J J (2tt) j 2tt x 02 
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To get the last line, we have used the inverse Fourier transform to represent the wave function in the coordinate 
space $'°)(z,r). Since the expression with the modified Bessel function QKi(Qxq2) depends on k ± we do not get a 
delta function 5^ (r — a; 01 ) after the integration over k 1 . This means that the parent dipolc recoils or changes its size 
because of the emission of the daughter dipole. If we insist that the recoil is small we can still make an improvement 
over the original LL dipole formula by making the following approximations 

(72) 

and after performing the k\ integration which in this approximation gives delta function we finally get 

2gt a ^°\z,x 01 ) ±-Q 0lKl (Q 0l x m )^^ . (73) 

27T XQ2 

This is an amplitude for the one gluon emission in the coordinate space improved by taking into account the next, 
subleading term in the energy denominator. Note, that it is a simplified and special version of (|29[) . The difference 
is that we do not keep all the terms in the energy denominators but only the ones which are related to the 'parent' 
emitter and the last gluon in the emission. This allows to factorize the splitting from the rest of the cascade. Therefore 
in this case each emission will be governed by the same kernel (|73[) whereas in the exact case the splitting is governed 
by more complicated expression with the Bessel function whose order changes with the number of the gluons in the 
wave function. 

The last expression (|73[) obviously reduces to the original LL dipole formula, compare (|59[) , by expanding the Bessel 
function K\ for the small values of the argument 

Qoi^i (Qoi ^02) = — Ki (—Vz) * —, for —VI -> . (74) 
The expression (|73|) becomes in this limit 



2.<$ (0) feoi^)^^, (75) 



which is the original LL formula [l[ as expected. Therefore, ([73")) is an improvement over the original formula (|75|) . 
Note that, there is a similarity of (|T3"|) with the form of the dipole formula for the F 2 structure function 

F T , L (x,Q 2 ) = -f— f d 2 r [ dz \*l*f(r, z,Q 2 )\ 2 a(r,x), (76) 
47^ a cm J J 

where ^ q is the wave function for the splitting of the virtual photon into a qq pair (dipole), and a is the imaginary 
part of the forward scattering amplitude of the qq dipole on the proton, called the dipole cross section, which describes 
the interaction of the dipole with the proton. In addition, r is the transverse separation of the quarks in the qq pair, 
and z is the light-cone momentum fraction of the photon carried by the quark (or antiquark). As usual, — Q 2 is the 
photon virtuality and x is the Bjorken variable defined as x = Q 2 / s with s the total energy of the j*p system. The 
wave function of the virtual photon is given by the following equations: 



,99|2 



3 a 
~2 



^Lj^e 2 {[z 2 + (l-zm 2 f K 2 (Q f r)+ m 2 K 2 (Q f r)} , (77) 
/ 

l*T*l a = ^KW-W}. 
/ 

where the sum is performed over quarks with flavor /, charge e/, and 

Q) = z(l - z)Q 2 + in} . (78) 
We see that whenever z or 1 — z are small the wave-function can be approximated by 
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in a close analogy with the dipole kernel in the LL approximation, compare the square of (|75|) . We see that ([73)) 
contains the kinematical constraint in the transverse space which is realized by the exponential tail of the Bessel 
function for large values of its argument 

QviKi{Q m x 02 ) = —Ki[—^fz) ~ J exp(- — y/z] , — y/z -> oo , (79) 

xqx \x Q i J y 2XQ2XQI \ x i J x 01 

We have deliberately chosen the notation of Q variable in (p9")) to make contact with the analogous variable in the 

wave function of the photon Q = Q 2 z(l — z) (|77I78|) . Of course we still have here small z, so that (1 — z) ~ 1. 
The improved dipole kernel coming from squaring (|73[) plus the term from the second graph, reads 



Xr\0 TT rr /-TT 



£02 



- Qoi K l(Qoi x 12) Z 



X12 



d X_2 Qqi 



Kf(Q 01 x Q2 ) + Kf(Q 01 x 12 ) - 2K 1 {Q m x 2)K 1 (Q 01 x 12 ) -° 2 ' ~ 12 , (80) 

X02X12 



with 

Q Q1 = — yfz , (81) 

compare Eq. (|72[) . We will refer to it as the quasilocal case because in this approximation we are keeping only terms in 
the energy denominators which refer to the daughter and the parent dipole, without any other dipoles in the cascade. 
It is straightforward to verify that the (|80|) simplifies to (|6Tj) when z — > 0. The kernel (|80|) is very similar in form 
to the one with the massive gluon, which also can be expressed in terms of the Bessel functions. Here, however the 
argument of the Bessel functions depends on the Q 01 , and consequently on the longitudinal momentum z, see (|5Tj) . 
The transverse and longitudinal momenta are not separated any more, even though we can still use a single, closed 
integral equation for the evolution of the dipole amplitude in the rapidity. In the LL approximation the evolution 
depended only on the previous step in rapidity, with the branching that was independent of the rapidity or x. The 
modified kernel (|80p contains the branchings which depend explicitly on the longitudinal variable, and therefore on 
all the steps in the evolution in rapidity. This is a qualitative difference as this means that there is now a 'memory' 
in the evolution of the system of dipoles. The probability of the emission of next dipoles depends on the evolution 
variable ('time') z. 



C. Soft gluon emission with the hard gluon displacement 

The approximation that we made in deriving (|80[) preserved an important feature of the soft gluon limit: the 
variation of the emitter gluon position in the process of soft gluon emission is neglected. Thus, transverse position of 
the harder of the two daughter gluons coincides with the transverse position of the parent gluon. A more accurate 
treatment of the emission process shows, however, that it is not the case. In order to demonstrate is explicitly, we go 
back to expression (|71[) with Q = k x y/z, 

2gt« f d 2 r*W(r,z) [ ^e* ^^) ±.Q Kl (Qx 02 )^™ , (82) 
J J (2-n-y 2n x a2 

and perform the integral over Cpk^ to get, 

2gt a J d\& Q \r,z) J ^^-^o(fcikoi -rD^^^^Cfcx^os)^ 



XQ2 



2gt" I d\$M(r,z)^M±^ f ^Uk^-r^K^yfzx^) 



x 02 2tt J 2n 



J ~ Vj x%2 (2tt) 2 [{x m ~ry + {y/z^x m y 



(83) 



The last term on the r.h.s. is a Cauchy type distribution in two dimensions which gives delta function S^ 2 '(x 01 — r) 
when the scale parameter y/zxo2 — ► 0. In this limit one reproduces expression Q75|) . Last line of (|83p gives more a 
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precise expression in the situation when the change of the emitter gluon position is taken into account. Note, that 
the integration in (153)) over the transverse vector r occurs already at the amplitude level, so the evolution equation 
for the dipole scattering amplitude with this kernel would involve a cumbersome triple integration over transverse 
coordinates. Finally, let us conclude that in a general case of the gluon splitting, the transverse position of neither of 
the daughter gluons coincides with the position of the parent gluon. 



D. Modified dipole evolution kernel and the NLL BFKL kernel 

In the previous subsection we proposed an approximate scheme to go beyond the soft gluon limit in the description 
of the onium wave function. In this scheme one better accounts for the kincmatical effects in the evolution of the gluon 
wave function. A significant part of the phase space for the gluon emissions is cut out. The effect is of course subleading 
with respect to the leading logarithmic approximation or the original dipole cascade. Nevertheless, the constraint 
introduces important corrections and cannot be neglected in the analysis of the color dipole / BFKL evolution beyond 
the LL approximation. In particular, we shall demonstrate explicitly that in the NLL approximation the modified 
dipole kernel (|B"D|) saturates the double collinear logarithmic contributions that are found in the NLL BFKL kernel. 
As discussed in [4l[ these contributions violate the renormalization group constraints and should be resummed. 
The resummation of these spurious terms within the conventional i-channel formulation of the BFKL evolution in 
the momentum space is conveniently realized by imposing the so called kinematical constraint (or the consistency 
constraint) [H, [33, 53] ■ The modified dipole kernel ([50]) and the BFKL kernel with the kinematic constraint have 
similar origin and describe the same effect, a modification of the evolution kernel due to a more accurate treatment 
of the kinematical effects. In what follows, we shall discuss in detail the problem of the double collinear logarithms 
in the dipole representation. 

In the exact next-to- leading calculation [l(J, [ll|, [l2|, [IH, [H| the terms with the largest power of the collinear 
logarithms come from the following part of the NLL BFKL kernel in the coordinate representation, 

X™f®N Y = -2al f ^ -£°V log ( ^) log ( [Ny(x , x 2 ) + N Y {x^x 2 ) - Ny^x,)} , (84) 
J 2tt Xq\ 2 x\ 2 \xoiJ \xoiJ 

where Ny{x i ,Xj) in the above equation is the scattering amplitude of the color dipole (j£i,Xj) with the target and 
Y = ]nl/x is the rapidity between the incoming dipole and the target particle. It is important to note that these 
double logarithmic terms are the same in the QCD as in the conformally invariant N = 4 SYM theory [TJ. Thus, we 
conjecture that this universality originates from universal kincmatical effects in gluon emission. 



FIG. 3: Dipole configurations in two collinear limits. Left plot: X02 <S xqi and X12 ~ XoV, right plot: X02 3> £01 and x\i 
X02 ■ 



Let us compare the behavior of (|80[) and (|84[1 in the collinear and the anti-collinear limits. These are limits when 
the parent and daughter dipole sizes are strongly ordered. In the anti-collinear limit one of the daughter dipoles is 
much smaller than the parent dipole and we have, 

(x 2 < X01 an d £12 ~ X01 ) or (^12 < X01 and x 2 ~ xqi ) . 

This situation is depicted on the left plot in Fig. [3j In these configurations no super-leading logarithmic enhancement 
of NLL BFKL kernel occurs, in particular the part of the NLL kernel |84|) is suppressed as 



log 



X12 
xqi 



log 



xqi 
x i 



or log 



X02 
xqi 



log 



xqi 
Xqi 







Thus, no anti-collinear double logarithms are generated. In the collinear limit both daughter dipoles are larger than 
the parent dipole. This configuration is shown in the right plot in Fig. [3l 



xq2 > xqi and x i2 ~ x 02 
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In this case the most singular part ([84]) of the next-to-lcading kernel tends to 

3CS£°«*r - log 2 (^) [-..], 

where [. . . ] stands for the expression in the square brackets as in ([51)) . 

Note that kernel ([84]) is large and negative in the regime where the produced dipoles x 02 and x 12 are very large 
compared to the dipolc xqi . This is exactly the same regime where the modified kernel (|80|) is exponentially cutting 
off the contributions of the large (02), (12) dipoles. Clearly, in such configuration the phenomenologically improved 
kernel ([50]) is better behaved than the NLL kernel. 

In the case of the kinematic constraint in the BFKL equation in momentum space, the double collincar logarithmic 
term is robust. These terms are entirely determined by the dependence of the effective cut-off scale for the transverse 
momentum on the longitudinal momentum fraction and does not depend on the details of the cut-off procedure. The 
same property should hold true in the position space. Thus, for the purpose of extraction of the contribution enhanced 
by the double collincar logarithms, it is sufficient to use a simplified dipolc kernel 

-^-9{x 2 1 -zx 2 2 )9{x 2 l -zx 2 2 ). (85) 

:C 02* Z '12 

This kernel introduces the cut-off on large dipole sizes consistent with the effective cut-offs given by the more accurate 
expression (|80p . but it is more convenient to analyze. The physical meaning of this kernel is simple: if the daughter 
dipolc size is small enough, the deviations from the LL kernel of the kinematic origin arc small, and if the dipolc 
sizes are too large, the kincmatical effects introduce a distortion, which leads to a suppression of the emission. This 
suppression, in particular, is visible in the Bessel improved kernel. 

The integral evolution equation for the dipole scattering amplitude with this kernel reads, 

N 01 (x)=N^+a a f —^l^k_e(zx 2 01 -xx 2 02 )d(zx 2 01 -xx 2 12 )[N 02 {z) + N 21 (z)~N 01 {z)}. (86) 

We shall go to the differential form w.r.t. the longitudinal variable x. One gets, 

o p j 2 2 ( 

-x—N 01 (x)=a s / -^^V 9{x 2 01 -x 2 02 ) 9(x 2 01 -x 2 12 )[N 02 (x)+N 21 (x)-N i(x)} 

OX J Z7T X Q2 X 12 I 



~2tt 



nl dz 

. z 



e{zxl 1 



XX 



12) +5 



L 02 , 



1 01 

c 12 , 



[zx 01 xx 02 ) 



x [N 02 (z) + N 21 (z) - N 01 (z)} 



(87) 



Clearly, due to limitations in the x 2 integration imposed by the 9 functions, the first line does not introduce the 
double logarithms. Let us focus on the part of the equation containing the integration over z. This gives, after the z 
integration 



d 2 x 2 



;E 02 X 12 



( x 02 x 0l)^( a '02 x 12) 



N 02 \x-¥ 

X, 



T 2 

x _02 

2 
' 01 



N 21 x^ 



Noilx^ 2 - 



+ 



(x 12 x 01 )9(x 12 x 02 ) 



' 01 



iV 02 [x -ii + N 21 [x -P _ JV 01 s " 



'01 



'01 



Being interested in the NLL accuracy we may expand the dipole scattering amplitudes with the shifted x argument, 
e.g. 



Nij (x Xq 2 I x\ x 



(89) 



using a logarithmic variable log(x), 



N tJ (xx 2 2 /x 2 01 ) ~ Na{x) + x^-N^x) log I ^ 



dx 



1 01 



(90) 
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The terms Nij(x) (without the differentiation w.r.t. x) may be combined with the first line of (|87|) to recover the 
LL dipole emission kernel. The emerging correction to the evolution equation takes the form, 



[8K®N] m {x) = g- / d 2 x 2 



'02^12 



8(Xq 2 x 0l)Q{ x 02 



f?a) log 



' T 2 ' 
x 02 



b 01 , 



x-^ [N 02 (x) + N 21 (x) -N 01 (x)} + 



(x\ 2 - x 2 01 )9(x 2 12 - x 2 02 ) log ( ^ ) x^- [N 02 (x) + N 21 (x) - N 01 {x) 



(91) 



Next, we shall use the BFKL equation to eliminate the x-differential, xd/dx Nij(x) , 



X dx~ Nij ( X>> = I d Xk 



2tt 



2 2 

X ik X kj 



[N ik (x) + N kj (x) - N i:j (x)} 



(92) 



Let us evaluate the action of the complete non-leading part of the expanded kernel on a probe function JVy (x) = 
i{x){x 2 a) 1 , for 7 — > 1. In this regime, the action of the LO order BFKL (dipole) kernel gives obviously, 



1-7 



N^(x). 



(93) 



This expression is inserted in (f9Tj) . First, we evaluate 



2n 



x 2 



'01 



^(^02-^01)^02-^12) lo S 



T 2 

x 02 
x 01 



1 



7 



N$(x) + N^(x)-N^(x) 



r(7) 



(7)/ 



(94) 



for 7 — > 1. The leading singularity at 7 = 1, is determined by the asymptotically large dipole sizes, 
and :rf 2 — ?~ 2 3> a^oi- Thus, the integral may be approximated by 



' 02 



r 2 » xl-y 



Us 

2n 



^ / ^^log| — 

c oi 



2a, 



,2-, 



"-oi 



(I-7) 3 



(95) 



where we took into account the effect of the constraint imposed by the function 



i)- 



Since x 12 



and \x 02 1 3> |iEoi 1 1 this 9 function constrains x 12 to a half- plane for which the angle between x Q2 and x 01 is smaller 
than ir/2. Therefore the angular part of the integration d 2 x 2 gives only ir and not 2ir. The same contribution as the 



one computed in (|95j) comes from the term proportional to 9{x\ 2 — x 2 )1 )9(3 



12 ' 



^02^ 



in (|91| . and the full answer reads, 



[5K®NW] 01 {x) = - 



^- 3 N$(x) 



(1-7) 







(1- 7 )S 



(96) 



This is precisely the same leading collinear singularity as the one found in the exact NLL kernel (|84|) , (Til. [T2I, [T3L [l4| . 

Finally, let us comment on the issue of the conformal invariance of the NLL BFKL kernel. The part of NLL kernel 
given by (|54"|) is scale invariant. This means that the rescaling transformation xj_ — > Ca^ leaves this part unchanged 
for an arbitrary number C. It is, however, straightforward to verify that is not conformally invariant in two 

transverse dimensions, for instance, the action of the kernel changes under the inversion x_ { — > =j. Interestingly 

enough, is the only part of the NLL BFKL kernel in N = 4 SYM theory that is non-conformal (of course in 

the QCD, already the scale invariance of the NLL BFKL/dipolc kernel is broken by the running coupling effects). As 
mentioned before the modified dipole kernel (|80p is also non-conformal in two dimensions since the arguments of the 
Bessel-McDonald function are exactly the same dipole size ratios — , — as in the above next-to-leading contribution 

J 1 X01 ' X01 ° 
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E. Diffusion in impact parameter space 



As is clear from the form of the modified kernel (|80p . the corrections from the energy denominators imply large 
modifications of the diffusion properties in the impact parameter space. The modified Bcsscl functions exponentially 
suppress the production of large size dipolcs above a z-dependent characteristic size 



2 x 



01 e -2x 02 Q 01 
02 



■VS 



2 ^02 X01 



r/01^02 



CO. 



The effective cut-off size ~ xoi/s/z grows with the decreasing z. This cutoff on the dipole size is analogous to the 
coherence effect in the cascade of the gluon emissions (64|. The maximal opening angle prevents the gluons from 
being emitted into a certain kinematic regime. Here, the effect is to prevent the emission of very large dipoles and as 
a result the diffusion in the impact parameter space is very much suppressed. 

To illustrate the diffusion better we inspect the Balitsky-Kovchcgov equation with the LL and the modified kernel 
respectively. We would like to stress however that the fact that we are using the modified dipole kernel in the non-linear 
equation is our assumption only, based on the analogy with the leading and next-to-lcading logarithmic calculations. 
This assumption is also motivated by the the fact that the dipole amplitude will still have the limiting value equal to 
unity in the saturation regime. In principle, the modification to the dipole kernel implies also the corrections to the 
triple Pomeron vertex, (tjq . l67j . This is an interesting problem which deserves further analysis but it is outside the 
scope of the present paper. 

Recall first the leading-logarithmic analysis from [68| where the power tails in impact parameter space were shown 
to emerge just after the first step in the evolution in rapidity of the BFKL and BK equations. We take the initial 
condition for the amplitude to be exponentially decreasing in impact parameter 



N(y = lnl/x = 0,x O i) = 2V ( 



where the impact parameter is defined as b 01 
with the LL kernel 



\ (a-0 + 3Li ) ■ 



(0) 
01 

Next, 



~ exp(-/i& i) , 

we perform one iteration of the dipolc-BK equation 



N 01 (x)=N^ + 



a s N c 
2tt 2 



dz 



d 2 x 



-2 _2 



[iVo2 + N 12 - Nqi - N 02 N- 



12 



(97) 



02' x '12 



For a small dipole xoi and at large value of the impact parameter boi the non- vanishing contribution to the amplitude 
comes from the configurations with large dipolcs, i.e. such that one end is at Xq and the other one is at x 2 located 
close to the center of the target where the field is strong. In that case 7V i ~ because of small dipole at large impact 
parameter and N02 ~ N12 ~ 1 and X02 ~ £12 ~ box- The r.h.s of the LL BK equation (|97|) is then 



d 2 x 



—2 ™2 



X 02 X 12 



[iV 2 + N 12 - A^oi - ^02^12] ~ -fi 



box 



(98) 



where 31 is the integration region where approximately X02 ~ box • Therefore in the case of the LL equation the 
diffusion in the transverse space leads always to the power like tails in impact parameter even with the exponentially 
falling initial conditions. This leads to the violation of the Froissart bound [f39l | even in the presence of the saturation 
corrections, as was first pointed out in (70l. l7ll. l72| . 

The situation changes when the modified kernel (|5D|) is considered. Here a similar analysis leads to 
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9{x m /\fzC - x 02 )- 



' 01 



+ &(x a2 - x i y/zC)- — ^ — exp 
2 x? 



[N 02 + N 12 - N 01 - iV 02 iV 12 ] , (99) 

b 02 \ x 01 /. 

where C ~ 1, and where we made approximations in the second term that xqi <C a;x2 ~ x 02 ~ boi which gives the 
dominant contribution. The first term is the short range contribution and the second one is the long range one. It is 
evident that the behavior of the scattering amplitude for small dipoles loi and at large values of the impact parameter 
boi is governed by 

. b i 



N(x i,b 01 ) — exp -2 ^fz 



(100) 
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Therefore there is an exponential but with the effective mass which becomes smaller as y/z/xoi decreases. One can 
actually see two limits in this behavior. We can take froij^oi fixed and so by decreasing z, or increasing the energy, 
the exponential tails become power like. One can also take 601 and z fixed, and change the dipole size xq\ . For small 
dipole sizes xm <^ &01 the tails are exponential, but for larger ones they become again power like. It is interesting 
to note that the largest effect of the modification is for the dipoles with small sizes, even though the cutoff inside 
the integral equation is acting on the large dipole sizes. This is result of the fact that the relevant parameter is 
XQ2/xoiy/z, i.e. is proportional to the ratio of the dipole sizes. The only way to eliminate the power-like tails is to 
put in (essentially by hand) the fixed mass term, which limits the range of the interactions (73j , for example 

Qoi ->■ Qoi + m2 ■ 

In this way one will get an amplitude which has exponential tails with fixed radius and this will lead to a behavior 
consistent with the Froissart bound [69| (modulo normalization). 

V. MHV SCATTERING AMPLITUDES FROM THE LIGHT-CONE WAVE FUNCTION 

A. Parke- Taylor amplitudes 

In the previous sections we have analyzed the (tree level) wave function with arbitrary number of gluons in the light 
cone formulation. This corresponds to the initial state evolution, that precedes the scattering. In this section we take 
into account also the final state evolution and evaluate the production amplitude of n-gluons in 2 gluon scattering. 
The exact tree level amplitudes with an arbitrary number of the external on-shcll gluons are known. These arc 
Parke- Taylor amplitudes Q (see Q for a comprehensive review) and can be recast in the following form 

M n = tr(t a H a2 ...t an )m(pi,ev,P2,e 2 ;...;p n ,e n ) , (101) 

{l,...,n} 

where a\, a2, ■ ■ ■ , a n , Pi,P2, ■ ■ ■ ,Pn and ei, 62, . . . , e„ are the color indices, momenta and the helicities of the external 
n gluons. Matrices t a are in the fundamental representation of the color SU(N C ) group. The sum in (|101[) is over 
the (n — 1)! non-cyclic permutations of the set {0, 1, . . . , n}. The kinematical parts of the amplitude denoted by 




FIG. 4: The n-gluon production amplitude with ( — , —,+,...,+) helicity configuration. All the gluons are taken to be outgoing. 

m(l, 2, . . . , n) = m(pi, e\;p2, €2] ■ ■ ■ YPm € n ) are color independent and gauge invariant. These objects have a number 
of important properties 

1. m(l, 2, ... ,71) are invariant under cyclic permutations of the set {0, 1, . . . , n}. 

2. m(n, n — 1, . . . , 1) = (—1)" rn(l, 2, . . . , n) (reversal symmetry). 

3. to(1,2, ...,n) + m(2, 1, ... ,n) + m(2,3, 1, . . . , n) + . . . + m(2,3, . . . ,n, 1) = (dual Ward identity). 

4. Incoherence to leading number of colors 

]T |M„| 2 = JV C "- 2 (JV C 2 -1) J2 {Ml,2,...,n)| 2 + 0(7V C 2 )} . 

colors {l,...,n} 
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It can be shown that the amplitudes where all the gluons have the same helicities, or only one is different from the 
others are vanishing (we assume that all the gluons are outgoing) 

m(±,±,...,±) = m( T ,±,±,...,±) =0. 

The non-vanishing amplitude for the configuration (— , —,+,...,+) depicted in Fig. 0] at the tree level is given by the 
formula 

m(l-,2-,3+,...,n+) = zg^ 2 — - — ( f } — -— — , (102) 

(12)(23) . . . (n — 2 n— l)\n— 1 n)(nl) 



where the spinor products are defined by Eqs. (jl7|) . (|18p . 

The comparison of the BFKL amplitudes with the Parke- Taylor expressions was also performed some time ago in 
[7-il |75| but only in the limit of the multi-Regge kinematics. Here, we shall re-derive the Parke- Taylor amplitudes 
within the color dipole picture in a much less restrictive kinematical limit. In order to compare the amplitudes 
obtained within the light cone perturbation theory formalism with the MHV amplitudes we need to consider the 
scattering process of the evolved wave function onto the target. We will simplify the problem by analyzing the case 
where the evolved projectile gluon scatters on a single target gluon which is separated from the virtual gluons in the 
projectile by a large rapidity interval. The exchange between the projectile and the target will be treated in the high 
energy limit. In this limit the interaction between the projectile and the target is mediated by an instantaneous part 
of the gluon propagator in the light cone gauge defined as 



DC" 



rfrf 



see for example [f|. We restrict the kinematics of the exchange, but still, the internal structure of the projectile gluon 
field is accurately represented. In principle, the technique applied could be used to evaluate the scattering amplitudes 
without any kinematical restrictions, but an analysis of the completely general case would be much more complicated. 




FIG. 5: Interaction of the wave function with 2 gluons with the target. The cross on the gluon line indicates that it is an 
instantaneous interaction. Vertical dashed lines denote the energy denominators; short dashes (blue lines), denominators in 
the initial radiation states; long dashed (red lines), denominators in the final state radiation. 



In order to demonstrate the equivalence with the MHV amplitudes it is essential to include the final state radiation 
in the wave function, which occurs after the interaction with the target has taken place. For example, in the case of 
the 2 — >■ 3 amplitude, we have to take into account graphs shown in Fig. [5J The gluon with the cross and labeled by 
the momentum I is the Coulomb gluon, and the gluon with momentum P is the target gluon. The latter one has a 
large P~ component. The dashed lines denote the energy denominators, which have to be taken before and after the 
interaction occurred. Similarly for the 2^4 amplitude we will need to take into account graphs shown in Fig. O In 
general for this kinematics the number of possible classes of diagrams is n — 3 with n — 4 energy denominators both 
in the initial and final state. Obviously within each of these classes one needs to sum over the different possibilities 
of the evolution of the wave function. For example in the case of graphs depicted in Fig. [5] there is a corresponding- 
set of graphs where the last emitted gluon is coming from the splitting of the most upper gluon in the hadron wave 
function. 

The dominance of the instantaneous (Coulomb) gluon exchange follows from the fact, that the other possible types 
of exchanges are suppressed w.r.t. the dominant exchange by large energy denominators. One can show easily that 
the relative suppression factor decrease with the total energy like M 2 /s, where M 2 is the invariant mass of the system 
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FIG. 6: Interaction of the wave function with 3 gluons with the target. The cross on the gluon line indicates that it is an 
instantaneous interaction. Vertical dashed lines denote the energy denominators; short dashes (blue lines), denominators in 
the initial radiation states; long dashed (red lines), denominators in the final state radiation. 



produced by the projectile fragmentation. At fixed M 2 , and for s — ► oo, the relative contributions of the non-Coulomb 
exchanges can be made arbitrarily small. Also, the contribution from scattering mediated by the four-gluon vertex is 
suppressed. 



B. Factorization of fragmentation tree amplitudes 



In constructing the full MHV amplitudes in the high energy limit we shall use the results for the initial state 
evolution supplemented by the analysis of fragmentation of the scattered states into on-shell final state gluons. 

The key feature of the fragmentation amplitudes is the independence of fragmentation of the gluon trees originating 
from different parents. We consider a (typically disconnected) fragmentation tree with a topology © — Oi U . . . U m 
where ©i,...0 m arc topologies of the fragmentation trees of m parent virtual gluons into n on-shell gluons, so 
that the first parent gluon, denoted by (12... ni) fragments into gluons (1, 2, . . . , n\), the second parent gluon, 
(ni + 1, ?ii+2, . . . , rix) fragments into gluons (rii+1, rix+2, . . . , n 2 ), and so on. Then the amplitude of the fragmentation 
tree 0, denoted by T& may be factorized into fragmentation amplitudes of the parent gluons, Te ; i n the following 
way, T e [(l • • • Hi), (ni + 1 . . .n 2 ), ■ ■ ■ , (n m _i . . .n) -> 1,2, ...n], 

Te[(l...ni),(ni + l...n 2 ),...,(n TO _i...n) -> 1,2, ...n] = T ei [(l...ni) -> 1,2, ...m] 

x Ts 2 [(m +I...TI2) -> ni + l,m + 2, ...n 2 ] x ... x T e ,„[(n m _i + 1 . . .n) -> n ro _i,n m _ 1 + 1, ...n], (103) 

provided that we sum over all possible light-cone time orderings of the splittings within all connected trees, while 
preserving the topologies. This property is not obvious in the light-cone formulation, as for a given time ordering of the 
splittings, the variables related to different trees are mixed in the energy denominators, and the factorization property 
(|103[) does not hold. It only holds if the complete summation over all possible splitting orderings is performed. In fact, 
the factorization property of independent tree amplitudes is quite intuitive. Nevertheless wc shall provide a formal 
proof of this property within the light-cone field theory framework. 

In order to prove (|103[) we shall consider two sequences of splittings, (ai), i = 0, 1, . . . , m and (bj), j = 0, 1, . . . , n , 
that occur in two topologically disconnected trees, a and 0b, after the scattering. The internal topology of each 
tree is arbitrary, in particular each tree, a and ©b, is not necessarily connected. The light-cone times of splittings et^ 
and bj will be denoted by r(ai) and r(bj) respectively. In the computation of the fragmentation amplitude of the tree 
= a U ©b one sums over all possible time-orderings of splittings. We shall denote the splitting ordering by T(O). 
This summation may be decomposed into three nested summations: two independent external summations of the 
time orderings within each tree, a and ©b, allowed by the topology, and the internal summation over all orderings 
of splittings between trees a and 0b, with a fixed ordering within each tree, 

E = E E E • 

T(e a ue b ) T(@ a ) T(e b ) T(e a ue b ) | T(e„)AT(e 6 ) 

The two external summations are performed for each tree independently, so they preserve the factorization property 
(11031) . Therefore, it is enough to prove Q103P for the arbitrary fixed splitting time orderings, T(0 a ) and T(0b). Thus, 
in what follows we fix the topologies a and ©b and the splitting orderings T(0 a ) and T(©b). It is convenient to label 
the splittings according to the time order, r(ai) < r(a2) < . . . < T(a rn ) and t(&i) < t(&2) < •• • < t(6„). 
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When the topology of the tree is fixed, then the contributions of splitting vertices to the amplitude docs not 
depend on the time ordering of the splittings. Furthermore, this contributions factorize into a product of independent 
contribution from trees a and Ob- This is not true only for the energy denominators that, in general, mix the variables 
of both trees. So, in order to complete the proof it is enough to show that the sum over relative splitting orderings, 
T(0 a U Ob) | 7(0 a ) A T(8b) of the products of the energy denominators of tree factorizes between the trees 9 a 
and 0f,. The energy denominators for tree a , in the absence of tree 0b, will be denoted by A\, A2, . . . , A m , where, 
in the final state we choose to attribute the denominator to the line before splitting. Analogously, the denominators 
for tree Ob, in the absence of tree O a , will be denoted by B\, B-x, ■ . - , B n . We shall prove the following factorization 
property for energy denominators: 



E 

(i p ,p— l,...,n+m) 



n+m 

n 



A, 



1 



Br, 



m n -. 



(105) 



where we used the additivity of the energy denominator of two evolving trees and introduced auxiliary symbols 
Aq = Bq = 0. The series (i p ,p — 1, . . . ,n + m) represents the ordering of splittings between the chains in the 
following way: i p is the number of splittings of a corresponding to p splittings of a U 0& assuming the time ordering 
T(6 a U 9b) I 7(0 a ) A T(0f,) , and p — i p is the corresponding number of splittings of Ob- 

Equation (|105[) may be proven by double induction, with respect to m and n. Let us first take n = m = 1. Then 
we have: 



in accord with (|105|) . Next, for 
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1 



1 
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Ai A 1 + B x Bi A x + Bt A X B X ' 
1 and an arbitrary m > 2 wc consider the l.h.s. of (|105[) : 



(106) 
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Ax +B-L A q + Bx A q A m Bi Ai + Bi A m + Bi 




Ax + Si V A 2 A m Bi Ax Am-xAr, 



, (107) 



where the first two terms in the square bracket in the middle line of (|107[) were evaluated using (|105[) for m — 1. This 
proves (|105[) for n = 1. 

In the second step we use the mathematical induction w.r.t. n, for a fixed m. The first step is to split the sum, 

E 



into contributions Si, S2, . . . S m , S m +i in which the first splitting in the 6-sequence, 61, occurs 
before 01, between ai and 02, . . . , between a m _i and a m and after a m , respectively. In each of these contributions 
we sum over all the allowed orderings of the later splittings, 62, . . . ,b n , w.r.t. the at splittings. Clearly, when 61 
occurs between a q and a q+ i (1 < q < m — 1), then all the later splittings, 62, ■ ■ ■ , b n , must occur after the splitting 
a q 

sequence a±, 02, 



It implies that in E g it is necessary to perform the complete summation over 62, . . . , b n orderings relative to the 



The splittings {ai,a 2 , 



,a q ,bi,b : 



2, . . • , b n _i} (with all possible orderings) are then followed 
The summation of the products of energy denominators 
b n -i}, preserving the ordering of (a,) and (bj) in T, q may 
be evaluated by applying (|105[) for the number of a-splittings equal n — 1 and the number of 6-splittings equal q. The 



by the (ordered) sequence of splittings: (6 n ,a g +i, 
over all orderings of the splittings {ai, a 2 , . . . , a q , bi, 62, 



remaining denominators, corresponding to (b n ,a q+ i, 
respectively. The final answer for E g reads, 



,) read, Aq+B u A q+1 + Bi, 



A r 



+ B U A r 
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(108) 
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The l.h.s. of (| 1 05[) is equal to sum over contributions S g : 
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Ai + Bi ' 



(109) 



where the contributions resummed using (|105|) are given in the square brackets. All the terms in the r.h.s. of (11091) 
have a common prefactor lj(B 2 . . . B n ) that can be factored out. The summation of remaining factors in the l.h.s. 
of (|109[) is precisely the one evaluated in (|10T[) . and the result of this summation is given by l/(BiAi . . . A m ). This 
completes the proof of the induction step in n, with fixed m. Taking into account the obvious symmetry between m 
and n, it proves also the induction step in m, with fixed n. Therefore (|105[) holds for all m and n. This, in turn, 
proves (|103|) . 



C. Fragmentation of a single gluon 



For the comparison of the amplitudes obtained in the light cone perturbation theory with the known expressions 
for the MHV amplitudes it is necessary to compute both the initial and final state emissions. This is essential, as we 
need to take into a account graphs which include both type of processes as depicted in Figs. [5] and [51 We consider here 
the fragmentation of a single, off-shell gluon (labeled by (12 . . . n)) into the final state of n on-shell gluons, 1, 2, ... n. 
The n final state gluons have transverse momenta k 1} . . . , k n and the longitudinal momentum fractions Zi, . . . , z n . 
The initial gluon has transverse momentum n ^ and the longitudinal fraction Z(i... n ) where we again used shortcut 



notation k 



(l...n) 



J2j=ikj and 



'(12.. 



S"=i z i ( see Sec. 



We will denote the fragmentation part of the amplitude for 1 to n gluons as T[(12 . . . n) — > 1,2, ... ,n]. Let us 
consider the n = 2 case. Using the computational methods developed in Sec. IIIII we can write the expression for the 
fragmentation in the following way 



T[(12)-»l,2] 



9 



12 



fJ 



4 /2 M 2 



£-3/2 



1 

Vl2 



(110) 



To obtain the above expression we have used the formula for the 3-gluon vertex ([B]) and the definition of the complex 
variable vi 2 ([7]). The energy denominator for the fragmentation D 2 was obtained from expression (|12[) . The calculation 
for the fragmentation of one to three gluons reads 



r[(123)-l,2,3] = -^.g 2 

^3 



J l(23) 



t l/2 ,3/2 
?1(23)V23 U 23 



(12)3 



t l/1 ^3/2 
£(12)3^12 U 12 



1 



2 1 -(123) \ ^'1(23)^(23)^12 +^( 12)3 ^(12)3^23 
" / ^12^23 



ZlZ 2 Z 3/ 



where the energy denominator for this case 



D, 



E 



k} (fc(123)) 2 



Zi 



2(123) 



(112) 



It is straightforward to verify that the numerator in the last line of (| 1 is equal to the denominator D 3 . Note, that 
this is a case of n = 2 of the generalization of the formula Eq. P2"]) derived in Sec. IIIII in the kinematics when the 
transverse momentum of the first particle in the fragmentation chain is not vanishing fc(i23) ^ and when z^i 2 n ) !• 
We also used the fact that 



£l(23)&23 — £(12)3^12 



2(123) 
Z1Z2Z3 
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Based on these two examples for n — 2 and n = 3 we postulate a general formula, that generalizes computations 
carried out explicitly for n < 4. The conjectured amplitude of fragmentation depicted in Fig. [7] reads 



r[(12...n)^l,2,. 



2(12. 



ZlZ 2 



3/2 



1 



(113) 



It is interesting to note that the form of the above amplitude for the fragmentation is dual to that of the gluon wave 
function with n-components. Namely, upon the replacement of the gluon light-cone velocities Vij in (| 1 13[) by the 
gluon transverse positions , we obtain the formula (152p for the on-shell gluon wave function in the coordinate space 
(modulo the overall factors with z, fractions). 





FIG. 7: Pictorial representation of the fragmentation amplitude T[(12 . . . n) — * 1, 2, . . . , n] for the single off-shell initial gluon. 

In order to prove the above formula (|113p we will need first to show that the relation (|112p is valid for arbitrary 
values of n. The generalization of Eq. (j42|) can be easily shown by noting that in this case we have 



(i...i)(i+i...n)£(l...-t)(i+i...n) - E k 



Ft ' -(■•••») Ft 



/ jkj = - — J-«(l...n) > 



2(1. ..77) 



(114) 



which is also a generalization of relation (|4"3"1) for the case when 2j=i kj ^ 0. The analogous proof as in (|4"2")l can be 
shown to proceed as follows 



5j«ii+l^(l...t)(i+l...re+l) w (l...i)(i+l...re+l) — E E fcj 



71+ 1 



/^7 



,•=1 Z (l...n+D j=1 

71+1 



2(1. ..7) 



h 2 h 2 h n+1 



i — i J — J- 



^77 + 1 



1 k. 1 



^~i Zl ^ I~ ^ Zj ~ z n+1 4- 1 h ~ f-f ^+7 ^ 

7—1 7—1 7 — 1 7 — 1 7—1 7 — 1 



^7+1 



Again, the second and fourth terms in the [. . .] cancel and we are left with 

mj+1 



77+1 , 2 , 71+1 V^ n+± h 

Ekj_ _ k n+1 , _ 2~,j=i K 3 



t 2? Z n -]-l 

7=1 J = l 



E 71-|- J 
7 = 1 



J=l *J 



i=i 



^71 + 1 

2n+l 



n+1 

E z j ~ Zn+i 



2^7 = 1 2 



^j=l -J 7=1 ~' s 2 (1 . ,.„ +1) 



(115) 



y fc L _ (fc(l..,i+l)) = 



The proof of the conjectured amplitude (| 1 13(1 can now be performed by mathematical induction, using the fac- 
torization property (|103p and the relation (|116p . Let us take the fragmentation into n+1 gluons denoted by 
T[(l, 2, . . . , n + 1) — > 1, 2, . . . , n + 1] and represent it by lower fragmentation factors T[(l . . . i) — > 1, . . . , i] and 
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T[(i + 1 . . . n + 1) — > i + 1, . . . , n . + 1] by separating out the first splitting (in all possible realizations). We obtain 



T[(12... ra 2,.. .,n+l] 



E 



'(!.. A)(i+X...n+X) 



i)(»+l...n+l) 

x T[(l...*)^l,...,i]T[(i + l...n + l)->i + l > ...,n+l] I . (117) 



This expression is the final state analog of formula p5[) for the iteration of the wave function and it is schematically 
depicted in Fig.[5J This may be rewritten using formula (| 1 1 3(1 for lower number of gluons: 




FIG. 8: Pictorial representation of one step in the fragmentation, Eq. (|117[1 , 



n 

T[(12...n+l)-l,2 > ...,n + l] = g n ^ £ 



z (X...i)Z(i+X...n+X) 
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(Ul2 . . . Vi-Xi) (vi+Xi+2 ■ ■ ■ Vnn+X) 



(118) 



The expression 



z (X...i)Z(i+X...n+X) 



3/2 
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C 3 / 2 \ ■ ■ ■ Zn+X 

S(l...i)(i+1...„+1) 



?(12...n+l) 



3/2 



Zx ■ --Zi Zi+X ■ ■ ■ Z n +x , 

is independent of « and can be factored out in front of the sum. Taking the common denominator we finally arrive at 



r[(12...n + l)->l,2,...,n+l] = 



Z{X2. ..n+X) 
ZxZ 2 ■ ■ ■ z n+1 



3/2 
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D n +i 



Z)t=l C(l...i)(i+l...n+l) v *X...i)(i+X...n+X) V H+1 
VX2V23 ■ ■ ■ Vnn+X 



Z (X2. ..n+X) 



3/2 



ZxZ 2 ■ . -Z n +X J 1'X2"V23 ■ ■ -Vnn+1 



, (119) 



where to get the last line we used the proof in Eqs. (|115( II 16[) that the sum in the numerator in the first line is equal 
to the energy denominator D n +x- The above expression is identical to (|113|) for n + 1 which completes the proof. 



D. MHV amplitudes from the LCPT 



It is interesting to see how the MHV amplitudes emerge in the LCPT framework. We have derived the exact 
expressions for the light cone wave functions in Sec. IIIII and the fragmentation amplitudes of the virtual gluon in 
Sec. IV CI In this section we compute the 2 — > n + 1 gluon scattering amplitude in the limit when a projectile gluon 
produces n gluons off a single gluon target, and the produced states is far in rapidity from the target gluon. The 
situation is depicted schematically in Fig. The result will be shown to reproduce the known analytic form of the 
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FIG. 9: The 2 — » n + 1 on-shell gluon amplitude in the high energy limit. The gluon dissociates into the gluon cascade 
(indicated by a blob) which interacts via high energy gluon (with a cross) with the gluon a — > 6. The large rapidity difference 
Ali ~ AY2 between the light cone cascade and the lower gluon is taken. The arrows indicate the momentum flow: the gluons 
0, a are incoming and 1, . . . ,n,b are outgoing. All the gluons have + helicity and it is conserved. 



MHV amplitudes in the adopted kinematic limit. In fact, due to this kinematic limitation we do not provide here 
the exact derivation of the MHV amplitudes in the general case. Note, however, that the high energy approximation 
is made only for the coupling of the incoming and outgoing target gluon, and no other approximations are made. 
Therefore, the comparison of our results to the MHV amplitudes provides a thorough cross check of our approach and 
the calculations. 

We label the gluons in the way shown in Fig. [5J The incoming gluon labeled 0, develops into a cascade of gluons 
1, . . . , n and scatters on a target gluon a — ► b. For the high energy case the dominant contributions are given by the 
instantaneous exchange of the Coulomb gluon. The exchanged gluon can be attached to any gluon in the cascade. 
Therefore we include both initial and final state emissions. However, since we are considering color ordered amplitudes 
it will be sufficient to take into account the attachment to the gluons in the form of the diagrams shown in Figs. [2] 
and [5] The gluons and a are incoming with hclicitics +, and all gluons through the whole cascade down to the final 
state, carry the positive helicity. 

In the simplest case when the upper part of the diagram is just single gluon — + 1 we have only 2 — ► 2 scattering 
and the helicity is conserved. The MHV amplitude for 2^2 scattering (| 102[) in this case 



\M(a,0 -> 



(aOY 



(a0) (01) (16) (6a) 



5 2 ^, (120) 



1*1 



where we have used the fact that the spinor products can be regarded as complex square roots of the Mandclstam 
invariants, sec formulae (|19[) . and in this case they are equal to 



1(^)1 = 1(16)1 = ^, I (01) I = I (ba) I = Vl*l • 

The explicit computation of the Coulomb gluon exchange in the LCPT in the high energy limit gives 

|M(o,0- 1,6)| = f ±, (121) 
^/z a z b z zi \t\ 

Obviously results (|120p and (|121|) agree up to the conventional factor 1 / J z a ZbZQZ\, used in the LCPT amplitudes 

In what follows, we shall factorize the amplitude of the gluon exchange, ~ g 2 s/\t\ , from the amplitude of the 
projectile gluon evolution down to the n-gluon final state, Thus, we write: 

M{a, 6, 1, 2, . . . n) = ~^9 2 ^ *n(l, % ■ ■ ■ n), (122) 

V z a z b I C I 

where the prefactor l/y/z a Zb is, again, a result of the convention adopted in the LCPT, and as usual the momentum 
conservation 5-functions represented by A^ n+1 ^ (recall that this factor is implicit in the wave functions) are factored 
out from the amplitude. In general, we have, 

M(a, -> 6, 1, 2, . . . n) = ■^/z a z b z z 1 ...z n M(a, -> 6, 1, 2, . . . n) (123) 
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Let us now consider the 2 — » 3 process depicted in Fig. [2] Using the light cone framework we can rewrite the upper 
part of both graphs in the following way 



§ 2 (1,2) = *i(12')T[(12) -> 1,2] + * 2 (1,2'). 



(124) 



where the first term on the right hand side of the above equation (fragmentation in the final state) corresponds to 
the upper portion of left graph in Fig. O and the second term (the initial wave function) to the analogous part in 
the graph on the right hand side. To be more precise we represent Eq. (|124| graphically in Fig. [TO] Note, that the 
prefactor z 2 , related to the Coulomb gluon coupling to gluons 2 and 2', cancels against the factors of X/y/zi, coming 
from the vertices connected by gluon 2'. The notation 2' indicates that we need to take into account the transverse 
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FIG. 10: Pictorial representation of Eq. 
exchange of the t channel gluon. 



124|l . The down-pointing arrows indicate the momentum transfer caused by the 
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FIG. 11: Diagrams contributing to the 2 
by the exchange of the t channel gluon. 



4 amplitude. The down-pointing arrows indicate the momentum transfer caused 



momentum carried by the exchanged gluon (the longitudinal fraction of the gluon momentum is neglected). It is 
equal to I = — k(i2)- For the process depicted in the last graph in Fig.[l0]we have fc 2 ' = fc 2 + fc(i 2 ) = —k\. Putting in 
the explicit expressions for the wave function 4' 2 (1, 2') Eq. (|38[) and T[(12) — > 1, 2] Eq. Q113|) we obtain 
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For the case of the 2 — > 4 process we have the following contributions 

$3(1, 2, 3) = *i(123')T[(123) -> 1, 2, 3] + * 2 (1, (23)') T[(23) -> 2, 3] + * 2 (12, 3') T[(12) -»• 1, 2] + * 3 (1, 2, 3') 

The first term on the right hand side corresponds to the sum of the diagrams A and B in Fig. 111! the second and 
the third terms correspond to the diagrams C and D, and the last term is the sum of diagrams E and F. Again, the 
explicit calculation using (f3"5| and (|113|) yields 
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(127) 
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In deriving (| 1 27[) we used the following relations 

1 ki 



= h 



zi(z 2 + z 3 ) £ 1(23 ) 



i v (12)3' = ^(12) 



1 _ ^(12) 

2 3 (^1 + Z 2 ) £(12)3 



as well as £1(23)^23 = £(12)3^12 = Z1Z2Z3. The case for 2-^5 proceeds in analogy. The corresponding expression for 
^4 can be found by taking the following sum 

* 4 = #4(1,2,3, 4')+ 

* 3 (1, 2, (34)')T[(34) ^ 3,4] + * 3 (1, 23, 4')T[(23) ^2,3] + * 3 (12, 3, 4')T[(12) -> 1, 2] + 
* 2 (1, (234)')T[(234) -► 2, 3, 4] + * a (l, (234)')T[(234) -> 2, 3, 4] + 

# 2 ((12), (34)')T[(12) -> 1,2]T[(34) -»• 3,4] + *i(1234')T[(1234) 1,2,3,4] , (128) 

where we have considered all possible attachments of the Coulomb gluon. The explicit calculations are straightforward 
but lengthy, with the result being 
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The expression for the general n can be found from the following formula 



(129) 




FIG. 12: Schematic representation of formula (|130[) . The down-pointing arrow indicates the momentum transfer due to the 
exchange of the Coulomb gluon. The part of the graph to the left of this arrow (before the scattering) is the wave function ^f m 
and the part of the graph to the right are the fragmentation amplitudes T, see Eq. (|130p . 

n 

#„(l,2,...,n) = # m ((l...ni)(ni + l...n 2 )...(n m _i + l...n)) 

m=l (l<rai<ri2<— <n m -i<n) 

x T[(l . . . m) — > 1, . . . , ni] T[(ni + 1 . . . n 2 ) — ► n\ + 1, . . . , n. 2 ] ... T[(n TO _i + 1 . . . n) — > n TO _i + 1, . . . , n] . (130) 

This formula, depicted in Fig. [T2] is obtained by taking into account all possible attachments of the exchanged gluon in 
the cascade. The additional assumption is the factorization of fragmentation of virtual gluon (|103[) . For the complete 
proof one has to use the explicit expressions for the fragmentation amplitude T given by (|1 13[) and the initial state 
wave function given by (|38[) . The general proof is quite lengthy therefore we give it in the appendix A. The final 
result is the following general form for the >F n for an arbitrary number of emitted gluons 

*„(l,2,...n) = . 9 - 1 ^^^= • (131) 

kx/Zx ^ZxZ 2 ...Z n Z\Z 2 . ..Z n V12V23 ■ ■ ■ Vn-Xn 
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FIG. 13: Pictorial representation of <? n . The down-pointing arrow indicates the momentum transfer caused by the exchange 
of the t channel gluon. The sum over all possible attachments to the gluons in the cascade is performed. 



We see that it is a generalization of the case for n = 2, 3, 4 given in (]125[) . (|127p and (|129[) . Note that, vqi = — ^ (as 
we have chosen the transverse momentum of particle to vanish k — ). This object is graphically represented in 
Fig. 1131 In order to facilitate the comparison with the MHV amplitude we use the relation (ii + 1) = y/ZiZi + \ m j+i 
to get 



* n (l,2,...n) = f 



1 



k 



(l...n) 



y/z lZ2 ...z n ^ (01) (12) (23)... (n- In) 



(132) 



We also need to include the Coulomb gluon exchange between gluon a — > b and the gluons in the cascade. Let us 
rewrite the s/\t\ part of the prefactor in (|122p . 



s 

V\ 



(aoy 



(aO) (nb)(ba)k 



(i...n) 



0(s°) 



(133) 



Collecting together (|122| . (|123|) . (|132| . and (|133j) . we obtain (up to a complex phase factor) the amplitude for 2 — > n+1 
scattering 



M(0;a 



,n;b) 



(aoy 



(aO) (01) (12) (n- In) (nb)(ba) ' 



(134) 



which is equivalent to the MHV amplitude (|102[) . The approximation sign "~" is used because the exchange between 
the gluon cascade and gluon a was computed in the high energy limit. 



Summary 

In this paper we have analyzed multi-gluon cascades keeping the exact representation of the gluon kinematics 
using the light cone perturbation theory. For the choice of gluon hclicitics corresponding to the MHV amplitudes, 
we considered the components of the incoming gluon light cone wave function with an arbitrary number of virtual 
gluons. We found a hierarchy of recurrence relation between the multi-gluon components that holds at the tree level. 
The hierarchy was solved exactly in the case of the on-shell incoming gluon. A compact exact form of the real gluon 
wave function is presented. Interestingly enough, the natural variables which appear in the computation are closely 
related to the spinor products which appear in the computation of the maximally helicity violating amplitudes. 

We also improved the color dipole evolution equation at small x by relaxing partially the soft gluon approximation. 
Thus, we included into analysis an entanglement between the longitudinal and transverse degrees of freedom, as it 
follows from the form of the energy denominators of the light cone perturbation theory. We obtained a modified 
kernel of the color dipole evolution in which the emission of large dipoles is exponentially suppressed above the size 
dependent on the longitudinal momentum fraction of the softer emission in the dipole / gluon splitting. The new 
kernel contains modified Bessel function of the second kind and it is no longer conformally invariant in two transverse 
dimensions. It also leads to a much slower diffusion in transverse space than the original LL dipole kernel. 

Furthermore, we showed that the modified dipole evolution kernel gives at the NLL accuracy the most singular 
pieces at the collinear limit of the exact NLL BFKL kernel, related to double logarithmic terms in the collinear limit. 
Notably, these are the terms of the NLL BFKL kernel which violate the conformal invariancc in the two transverse 
dimensions both in the case of the N = 4 SYM theory and in the case of QCD case. This supports a conjecture that 
these terms are related rather to universal kinematical effects than to dynamical details of the theory. 
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We have also investigated the scattering of the gluon, described by the light cone wave function on a target gluon 
and the final state emissions of gluons. We showed explicitly within LCPT that fragmentation amplitudes of different 
virtual gluons emerging from the scattering into real gluons are independent, despite an apparent entanglement. 
The fragmentation amplitude has been computed exactly of a single virtual gluon with the positive helicity into an 
arbitrary number of real gluons with positive helicities. Interestingly, we found a duality between the initial state 
wave function of the real gluon and the fragmentation amplitude of the virtual gluon. This result may hint on a 
deeper duality between the initial and final state evolution in the gluon cascade. 

Finally, we have constructed the 2 — > n + 1 gluon amplitudes in the MHV helicity configuration by scattering the 
wave function on the target gluon, followed by the fragmentation of the scattered state. We chose the kinematical 
situation in which the projectile gluon produces n gluons in arbitrarily large rapidity separation from the target gluon. 
In this limit we reproduced the Parke-Taylor MHV amplitudes. 
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Appendix A 

We provide here the proof of the explicit expression for the upper portion of the scattering amplitudes Q13ip which 
will be done via recursive induction. Let us begin with the formula (|130|) for the function ^> n introduced in Sec. [V] 

n 

#„(l,2,...,n) = ^ * m ((l...ni)(m + l...n 2 )...(n m -i + l...n)) 

m=l (l<ni<n 2 <...<n m -i<n) 

x T[(l . . . m) — > 1, . . . , n{\ T[(ni + 1 . . . n 2 ) — > ni + 1, . . . , n 2 ] ... T[{n m -i + 1 . . .n) — > n m -\ + 1, . . . ,n] . (135) 

This formula arises when we take into account all possible attachments of the exchanged gluon between the cascade 
of gluons and the target gluon. Therefore we need to combine the wave function (initial state) with the fragmentation 
function (final state radiation) in all possible ways. This is accounted for by performing the summations in (|135|) . 
The above sum was used to compute the explicit expressions for \P 2 , ^3, ^4 which are given by the Eqs. (|125ll27ll29p 
by using the formulae for the wave functions and the fragmentation functions. Here we prove that using the Eq. (|135p 
we get the following explicit result for $f n (|131[) in case of the arbitrary n number of gluons 

*„(l,2,...n) = ff «-ifc=> _]== 1 1 . (136) 

kl/Zl y/ZiZ 2 ...Z n ZiZ 2 ...Z n Vi 2 V 23 ■ ■ ■ Vn-ln 

As a first step we need to take into account the momentum transfer from the exchanged gluon in the explicit expressions 
for the wave functions. For example, the wave function ^„ with n = m (with the momentum transfer included) which 
appears in the sum (|135[) has the following form 

* B (l,2,...,n') = (-I)"" V- 



^ziz 2 ...z n £ (12 .. 

.n— X)n s(12...n— 2)(n— 1 n) ■■■ Sl(2...n) 

i 



v (12...n-l)n' w (l 2 ...ri- 2 )(n-l n)' ••• V l(2...n)' 



(137) 



where primes indicate the momenta of the wave function before the momentum transfer has happened as discussed 
in Sec. [V] and shown in Figs. [TOland llll From the momentum conservation (the transverse momentum of the initial 
particle is zero) we get following relations between the primed and unprimed momenta 



"W — — k(i...n-l)i fe(n-ln)' — ~ k(l...ra-2) ) ■■■ > &(2...n)' — ~ 
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Therefore in general we have 

fc (i-i) fc (j+l...n)' _ k (i... 3 ) 



V(l...j)(j+l...ny 



These relations are true for any wave function appearing in expression (|135|) . Using these relations we can rewrite the 
wave function (|137| as 

*„(l,2,...,n') = (-l)-V- 1 ^^ — ; — h , (138) 

yJZ\Zl ...Z n ki fc( 12 ) . . . K(i...„-i) 

which is equivalent to (|47p . For each term in the sum (|135[) we have the factor which depends on the longitudinal 
degrees of freedom in the following way 

/ \ 3/2 

/ z (l...n 1 ) z (n 1 + l...n 2 ) • ■ • z (n m _i+l...n) 



y/ z (l...ni) z (ni + l...n 2 ) ■ ■ ■ z (n m -i + l...n) \ Z\ . . . Z n 

1 / z (l...n 1 ) z (n 1 + l...n 2 ) ■ ■ ■ z (n m - 1 + l...n) 



Z\ . . . z n y Z\ . . . z 7] 



(139) 



Therefore we can factor out the term 1 / Jz"i . . . z n in front of the whole expression. This factor is purely conventional 
as it comes from our choice of the normalization of the momenta. The common factor l/z\Z2 ■ ■ ■ z n in the parenthesis 
() will be left under the sum for the time being. It will prove to be useful to use another notation for the momenta. 
We introduce auxiliary variables pu+i defined as follows 

Pl2 = ki, P23 = fc(12)j P34 = fc(123)) ■ ■ ■ i Pn-ln = ■ 

Using ()135[) , the new variables, explicit expressions for wave function (|138p and fragmentation amplitude we get for 
the general case n 



* n (l,2,...,7l) = 



g n 1 



yJZxZi ...Z n y Zx . . . Z n V 12 . . . V n ^i n 
z (23...n) 1 Z(12)Z(3...n) 1 z (l...n— 1) 



z 2- ■ - Z n pi 2 1>23 • ■ ■ «n-ln z l z 2 ■ ■ ■ Z n Wl2P23«34 ■ ■ ■ «n-ln ^1^2 ■ ■ ■ 2W-1 ^12 ■ ■ ■ V n -2n-lPn-ln 

z (l...i 1 ) z {i 1 + l...i 2 ) z (i-2 + l...n) 1 



Z 1 ...Z n V12 ■ ■ .Ptiti+1 ■ ■ -Pi2i9+1 ■ ■ ■ Vn-lr 



E 

^(I...ii)^(ai+l...i2) z (»2 + l...i3) z (»3 + l---") 1 



E 



(-iy 



z (12) 1 2(23) 1 z (n-ln) 1 



^l2 2 W12P23 ■ ■ - Pn-ln Z 2 Z 3 Pl2«23P34 • • -Pn-ln z n-l z n Pl2 ■ ■ ■ Pn-2n-lV n -lr 



i-ir \ 

P12P23 ■ ■ -Pn-ln 



where means that we need to take higher order sums with 1 < i\ < i% ■ ■ ■ < i m < n — 1. We see that in the 

expression above, each term has a form of 

z (l...i 1 ) z (i 1 +l...i 2 ) ■ ■ ■ z (i m + l...n) 1 



Zl . . ■ Z n V12 ■ ■ ■ ■ -Pi 2 i 2 + 1 Pi m i m + 1 ■ ■ ■ V n -ln 

with 1 < ii < 12 ■ ■ ■ < i m ^ ti — 1, where m can range from to n — 1. In other words we can regard (|140j) as a sum 
over the different possible insertions of two kinds of elements pi k i k+ i and f^^+i into the chain with n — 1 sites. We 
write down the above formula in the form with the common denominator, which is equal to 

T> n = z l z 2 ■ ■ ■ z n V12V23 ■ ■ ■ V n -lnPl2P23 ■ ■ - Pn-ln ■ (141) 
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The numerator is then equal to (we omitted the g n factor here) 

N„ = Z(l... ra )Pl2P23 • • -Pn-ln ~ 

[ z l z (2...n)Vl 2 P23 ■ ■ - Pn-ln + Z(12)Z(3...n)Pl2V23P34 ■ ■ - Pn-ln H h z (l...n-l) z nPl2 ■ ■ ■ Pn-2n-l^n-ln] + 

22 z {l...i 1 )Z{i 1 + l...i 2 ) z {i 2 + l---n) Pl2 ■ ■ ■ V ilil + i . . . V l2l2 + \ . . -Pn-ln ~ 



l<~i\ <i 2 <n— 1 

E z (l...i 1 ) z (i 1 +l...i 2 ) Z (i2+l— i 3 ) Z {i3+l— n) 
Pl2 ■ ■ ■ Vi 1 i 1 + 1 ■ ■ ■ Vi 2 i 2 + 1 . ■ ■ Vi 3 i 3 + i . . -Pn-lr 

i ^ , . , , ^ , z l • • • z n 

!<<ii<i 2 <i3<n—l 



12) z 3 • • • z nPl2^23 • • • ^n—ln ~\~ z l z (23) • • • ^n^l2^23^34 • • ■ Vn—ln ~t~ ■ ■ ■ ~t~ Zl • • ■ z {n—ln)^12 • • • v n—2n—lPn—ln\ 

+ (-1)" ZiZ 2 ■ ■ ■ Z n Vi2 . . -Vn-ln ■ (142) 

We assume that the formula (|136p is true for n. In that case 

N n = Z 1 p 2 3P3i ■ ■ ■Pn-lnk(l... n ) ■ (143) 

Let us take the case of n + 1 which reads 

Xn+l = Zn. ..n+l)Pl2P23 ■ ■ ■ Pn-lnPnn+1 

[ziZ(2...n+l)Vl2P23 ■ ■ ■ Pnn+1 + 2(12) z (3. . .n+l)Pl2«23P34 ■ ■ ■ Pnn+1 H h z (l...n) z n+lPl2 ■ ■ ■ Pn-ln^nn+l] + 

2(i...i 1 )Z(i 1 _|_i...i 2 )2(i 2+ i...„_|_i) P12 ■ ■ • Uj 1 j 1+ i . . . V i2l2+ i . . . Pnn+1 

l<ii<i 2 <n 

+ + (-1)" +1 Z\Z2 ■ ■ ■ Z n Z n+ iVi2 ■ ■ ■ Vn-lnVnn+1 ■ (144) 

Obviously in every term in Ji n +i there can be either p nn +i or v nn +i present. We collect these terms and obtain 

^n+l — Pnn+1 ( ) + z n+lVnn+l [ z (l...n)Pl2 • • • Pn— In + 

z lZ(2...n) v 12P23 ■ ■ -Pn-ln + 2(i 2 )Z(3...ri)Pl2V23P34 ■ ■ -Pn-ln H h 2(l...n-l)2nPl2 ■ ■ •Pn-2n-lW„-l ra 

_ ^ 2(l...t 1 )«(i 1 +l...t 2 )% a +l...n) Pl2 ■ ■ • "iiii+1 • • ■ ^i 2 i 2 + l ■ ■ -Pn-ln + 

l<i± <i 2 <n—l 

+ + (-l) n+1 ZiZ 2 ...Z n Vl2... Vn-ln] ■ (145) 

The terms in the squared parenthesis [. . .] combine to — N„ of (|142[) . We make the inductive step and use (|143[) for 
this expression. Using p nn +i = &(i...n) the numerator N n+ i reads then 

N n +1 = Pnn+1 {z(l... n +l)Pl2P23 ■ ■ -Pn-ln ~ 

[ziZ(2...n+l)Vl2P23 ■ ■ -Pn-ln + ^(12) Z(3...n+l)Pl2«23P34 • • -Pn-ln H 1" z (l...n-l) z (nn+l)Pl2 ■ ■ ■ Pn-2n-l«n-ln] + 

z (l... ii) z (ii + l... i 2 ) z (i 2 + l ...n+1) Pl2 ■ ■ ■ "^iiii+l • • ■ u i 2 i2 + l ■ ■ ■ Pn-ln — 

l<zi<22<n— 1 



Zn+lVnn+lZlP23P34 ■ • • Pn—ln } , (146) 

where the last term comes from making the inductive step as described above and we have factored out p n n+i- Now 
we combine terms in the parenthesis {. . .} which contain either p n -in or v n —i n - We then observe that the coefficient 
in front of zi n n+i)Vn-in will combine to Nn-i and we can again use (|143|) for n — 1. Therefore the term proportional 
to u Tl _i„ in {. . .} reads 

— Z(nn+l)V n -lnZlP23P34 ■ ■ - Pn-ln ■ 

Therefore all the terms now contain p n —i n and it can be again factored out. We see that we need to apply this 
procedure recursively n — 1 times, and each time we pick up a term which will be proportional to ztj„, n +i) v jj+i 
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(where j decreases with each step). Making the inductive step we see that such term will contain Pjj+i which can be 
factored out. After performing all these steps the result reads then 

N n+ i = P23 • • - Pn-lnPnn+l \_ z (l...n+l)Pl2 — Z\ Z( 2 . . .n+1) v 12 — Z\Z(Z...n+l) u 23 — • • ■ — ^l^(nn+l) u n-ln — ZlZ n +lV n n+l\ , 

(147) 

where all Pjj+i terms with 2 < j < n have been factored out. We can use 

— z l z (k...n+l) v lk — z l z (k+l...n+l) v kk+l = ^ z l z k v lk — z l z (k+l...n+l){ v lk + Vkk+l) = — z \ z k^lk ~ z l z {k+l...n+l) v lk+l j 

(148) 

to rewrite the numerator as 

N„+l = P23 ■ ■ -Pn-lnPrm+l [z(l...n+l)Pl2 - ZlZ2"12 - z \ z ZV\z Zl^nVln - Zl^n+lUln+l] ■ (149) 

Using px2 = &i and z\ZjV\j = Zjk\ — z\kj we find that the terms in the [. . .] combine to Zi &(i... n _|_i) which gives 

N„+i = Z1P23 ■ ■ ■Pn-lnPnn+lk(l... n +l) , (150) 

and this completes the proof for ^>„+i- Let us finally note that in the relations (|147|) to (|150[) were also used in 
the explicit calculations done for n = 2,3,4. It is important to note that these relations hold also in the cases 
where Z(x... n +i) 7^ 0, which is important as in the proof we recursively move to lower n where the partial sums over 
longitudinal z and transverse momenta k can take different values. 
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